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Abstract 

We prove two asymptotic expansions of the generalized scattering phases. These phases are 
generalizations of the Birman-Krein spectral shift function associated to pairs of perturbations 
of the Laplacians of asymptotically hyperbolic manifolds. The first expansion, of 'heat type', 
holds for all 'long range' metric perturbations of the Laplacian, whereas the second one is 
shown under a non trapping condition. 

1 Introduction and results 
1.1 Introduction 

In this paper, we define and study some properties of the generalized scattering phases associated 
to a pair of self-adjoint elliptic differential operators (Pq,Pi) on an asymptotically hyperbolic 
manifold X of dimension n. The manifolds and operators that we consider are of the same type 
as those considered in or in many other related papers, especially [3U]. Typically, P can 

be the Laplacian of the quotient of the hyperbolic space HP by some discrete group of isometries, 
with infinite volume, and Pi a non compactly supported perturbation of Pq. The precise definition 
of the operators is given in subsection 1 1.21 Let us first explain the terminology. 
The generalized scattering phase of order q £ N is the distribution £ q defined by 



jf /(A)d£ q (A)= tr(/(^)-£±^/(P^ j, f£C ( 



(1.1) 



provided, of course, that the right hand side makes sense. Here P e — Pq + e(Pi — Pq) and the 
left hand side stands for (£ q , /} if (., .} is the duality between test functions and distributions. 
This equality actually defines the derivative of £ q and we choose for £ q the unique primitive of £ q 
vanishing near — oo, which is possible if Po and Pi are semi-bounded from below. 

The distribution £i is well known: it is the spectral shift function introduced by Birman 
and Krein |2l EH- We recall that £i is, in general, a measurable function which is defined if 
(Pi +i)~ N — (Pq + i)~ N £ Si for some N large enough. Here Si is the set of trace class operators. 
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The point of considering £ q for q > 2 is that we can relax this trace class condition and 
consider operators such that (Pi + i)~ N — (Pq + i)~ N £ S q , the Schatten class of order q (for 
instance the Hilbert-Schmidt class if q = 2). We furthermore emphasize that, if q > 2, £ q is a 
priori a distribution and showing that it is a (smooth or continuous or even measurable) function 
on the continuous spectrum is not trivial. The distributions £ q have been introduced by Koplienko 
|.'{2j for bounded or 1-dimensional Schrodinger operators and studied in the higher dimensional 
case by the author (see [Hl[7| for more details). We also quote a similar approach considered in 
|28j . In HUG], £q was called spectral distribution by analogy with spectral function, however the 
name scattering phase is natural as well since we have shown that 



The precise meaning of formula (|1.2fl . which is well known for q = 1 if (Pi — Po)(Po — z) _1 is 
trace class 0EB], is explained in 0. We just specify that det q coincides in many cases with the 
standard Fredholm determinant Det q defined for perturbations of identity by elements of S q (see 
[281 158) ) and that we need such an extension of Det q since (Pi — Po)(Po — z)~ l is not necessarily 
compact in general. We also refer to the recent paper [5] where similar determinants with q = 1 
are studied for hyperbolic surfaces. 

We now recall the definition of an asymptotically hyperbolic manifold. A complete non compact 
Riemannian manifold [X, G), without boundary, is asymptotically hyperbolic if it is isometric, 
outside a compact set, to 



Here Y is a connected compact manifold without boundary, and g(x) is a family of metrics on 
Y depending smoothly on x E [0, e _i? ). More precisely, if S 2 T*Y is the vector bundle of bilinear 
symmetric forms on Y and T(S 2 T*Y) is the space of its smooth sections, we assume that g E 
C°°([0, e~ R ), T(S 2 T*Y)) and of course that g(x)\T p YxT p Y is positive definite for each p E Y and 
x E [0, e~ R ). We furthermore equip Y with the metric <?(0). We could actually consider manifolds 
with finitely many such ends, i.e. with (A, G) isometric, outside a compact set, to a finite union of 
manifolds like l[1.3[l but we restrict our attention to the one end case for notational convenience. 
As explained in [23 ESj j such a manifold (X, G) can be obtained from a compact manifold with 
boundary (Z, G), with a boundary defining function x such that dZ = Y = {x = 0}, by setting 
X = {x > 0}, the interior of Z, equipped with the metric G/x 2 . These manifolds are also called 
conformally compact manifolds and the most basic example is the hyperbolic space H„. 

The results of this paper are two asymptotic formulas for £ q which are similar to the heat 
expansion and Weyl formula for the eigenvalues counting function on a compact manifold. More 
precisely, these asymptotics are of the same type as those obtained in Euclidean scattering by 
[HI EI 113 112 E31 E2 EH| for g = 1 and 00 for q > 2. See also [El El M EH and |2EJ in more 
geometric frameworks. The scattering phases are natural and basic tools of scattering theory in 
view of their relation with time delay |47| . relative scattering determinants |51l71l43| or resonances. 
More specifically, their asymptotic behavior is of interest for several reasons such as relative index 
theory gj, trace formulas [T3E3EHI] or Breit-Wigner formula [2"2l |4H I5| . 

The most popular scattering phase is Birman-Krem's function £i but its use leads to restrictions 
on the pairs of operators as already mentionned. For instance in 0, the authors are able to define 
the determinant deti for a pair of operators which are, up to a unitary transform, the Laplacians 
associated with two metrics Go and Gi as above for which gi(x) — go(x) — 0(x 2 ). This last 
condition implies that (Pq + i)~ N — (Pi +i)~ N is of trace class since they work in dimension 2, but 



e q (A) = Urn A arg det q ((P — A - i5)(P - A - iS)- 1 ) . 



(1.2) 



((P,+oo) xY,dr 2 + e 2r g(e- r )) . 



(1.3) 
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in higher dimension, the difference (Pq + i)~ N — (P\ +i)~ N is not trace class, in general, under the 
sole condition gx(x) — go(x) = 0(x 2 ). Our theorem 11.11 combined with the method of |7] proves 
directly the existence of det q ((Pi — z)(Pq — z)^ 1 ) for q > n for any n under the weaker condition 
that g x {x) - g (x) = 0{x). 

For the proof of theorem 11.21 we adapt the method that we used in 0, namely a refined 
analysis of Isozaki-Kitada's construction We recall the principle of this method, which has 
only been used in the Euclidean context so far, in subsection 12.31 We devote section to the 
relevant estimates on geodesies in the hyperbolic framework and the explicit construction is given 
in subsection l5.ll 

1.2 Notations 

For any C q function T £ of the real variable e £ [0, 1], scalar or vector valued, we set 
We will also use the notation 

1 f 1 d^ 1 

Note that, if T e is a primitive of T e , we have [%] q = {T e }q. 

We will have to consider distributions smooth with respect to a parameter. We shall say that a 
family of distributions u T E S'(R) is with respect to r E J, an interval of R, if for all / E <S(R) 
the function (u T ,f) is C J on J. In particular, d{u T and JjU T dT are defined by 

(dlu T ,f) = di(u T J), (Ju r dr,f\ = J(i*r,f) dr. 

Another useful distributional notation is the following. If H is a self-adjoint operator on a separable 
Hilbert space TL and if T is an operator acting on (some subspace) of TL such that the map 
/ i > tr (f(H)T), with / E 5(R), defines a distribution, that is if the trace is well defined and 
depends continuously on / E 5(R), then we shall write this distribution 




if E(p) = E(—oo, fi) is the spectral projection of H on (— oo, /it). 

We will also use extensively Schatten classes S q of real order q > 1. We simply recall that, 
by definition, a bounded operator A on a separable Hilbert space Ti. belongs to S q = S q (7Y) if 
\A\ q = {A* AY' 2 is trace class and that the norm ||.|| q on S q is defined by \\A\\\ = tx(\A\ q ). We 
will need the following Holder type estimates 

P£||q<Pll qi ||£|| qa , Q-^q^+ta 1 (1-5) 

for all A E S qi and B E S q2 . This estimate still holds if A (resp. B) is bounded with q\ = oo 
(resp. (72 = oo), which is consistent with the notation ||.||oo for the operator norm on Ti.. For a 
more general presentation of Schatten classes, we refer to and |53| . 

On the manifold X, we will mainly work near infinity. If y — (yi, • • • , y n —i) are coordinates on 
Y defined on U n -i C Y, such that y : U n -\ -t!lc R™ _1 is a diffeomorphism, then r, yi, ■ ■ ■ , y n -i 
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are coordinates on an subset U n C X diffeomorphic to (R, oo) x ft. We call (R, oo) x ft a chart 
at infinity. We can write Y as the finite union of open sets U n -i and thus we get an atlas of a 
neighborhood of infinity on X given by a finite number of charts at infinity. If necessary, we can 
choose each f2 to be convex, since we fix the coordinates once for all in each chart. We assume 
that, in each chart at infinity, the volume density giving the L 2 structure on X can be written 

dvol = a(r, y)\dr A dyi A • • ■ A dy n -i\ (1-6) 

with a smooth and bounded on (i?, oo) x Q. In particular, this implies that the pullback on X of 
any function u G L 2 ((R 7 oo) x Q,drdy), supported in the chart, belongs to L 2 (X) — L 2 (X,dvol). 

We now describe the operators Pq and P\. Let us start with an example. The expression of 
the Laplace Beltrami operator associated to 1(1. 3|) in a chart at infinity is 

A G = -d 2 - e- 2r A^ -(n-1- e -r 9 ^ et 9( x ) | 

V 2detg(x) 1 

where A 9 ^ is the expression of the (negative) Laplacian of Y associated to the metric g(x). 
Furthermore dvol(G) — e^ 1 ^ \J 'det g(e~ r , y)\dr A dyi A • • • A dy n -\\ is the volume form induced 
by G, hence 

e (n-lW2 AGe -(„-lW 2 = _q2 _ e -2r A9 (e- n + £ 2 f -r d x det 9 ( X ) _ 

2detg(x) 

with c„ — (n — l)/2, is selfadjoint w.r.t the density y / det g(e~ r , y)\dr A dyi A • • • A dy n ^\ \ which is 
of the form H1.6(l . Guided by 1)1.7(1 , we shall consider operators Pq , Pi which are both second order 
elliptic differential operators, symmetric w.r.t to dvol, whose expressions in each chart at infinity 
are 

Pi = -d 2 r -g J (e- r ,y,e- r d y )+e- r ]T vf (e^, y)^^)^ j = 0, 1. (1.8) 

\a\+l<l 

Here gj(x,y,ri) is the principal symbol of A g ^ x \ j — 0, 1, that is the expression of the metric gj on 
the fibers of T*Y with coordinates r)i, ■ ■ ■ , r\ n -\ dual of yi, ■ ■ ■ , y n -\, and the functions v"' l (x, y) 
are smooth and bounded on [0, e~ R ) x ft. The fact that Pi is a perturbation of Pq is reflected by 
the assumption that 

9i\x=o = 9o\x=o- (1-9) 
We will use extensively the principal symbol of Pq + e(Pi — Pq), denoted by p ei which has the form 

Pe(r, y, p, v) = P 2 + e~ 2r 5 £ (r, y, 77) 

where g e (r,y,r/) — go(e~ r ,y,i]) + e(gi — go)(e~ r , y, 77) and p is dual variable of r. Note that 
9i{x,y,V) -9o{x,y,t]) = C(a;)|r7| 2 by (JTU). 

1.3 Results 

In the next two theorems, Pq and Pi are two operators as described above such that (|1.9fl holds. 
We recall that P e = Pq + e(Pi - Pq). 
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Theorem 1.1. i) For all q > n and f G [/(P c )]q is trace class and there exists a unique 

£ q £ 5'(R) which vanishes below inf (cr(Po) U <r(Pi)) such that 

//(A) d£ q (A) = tr[/(P e )] 
Jr 

ii) The Laplace transform o/£ q has a complete asymptotic expansion as t j 0+, namely 

tr[e- tPe ] q -^ n/2 ^a fe t fc , with a =r(|+l) {2n)- n uj n j [dvoZ £ ] q . (1.10) 

k>0 ^ X 

Here dvol e is volume density obtained naturally from p f , that is in local coordinates 

dvol t = e^-^det (d*g e (r,y, r,)^)^' 2 \drAd Vl A---A dy n -i\- 
The other coefficients a\,a2,-- m can be expressed as integrals of functions of the symbols of P t . 

Example. If (X, Go) and (X, G\) are of the form (|1.3|) . outside a compact set, associated respec- 
tively to g and g\ satisfying 11.9fl . then the operators Pq — e (™- 1 W 2 A G °e~(™~ 1 ) r / 2 — c 2 and 
Pi = [/i/2 e («-i)r'/2 A G le -(n-i)r/2 [/ -i/2 _ c 2 } with jj = dvolid) / dvol(G a ) , satisfy the assump- 
tions of the theorem. 

The leading term of (|1.10|) involves the regularized volume term oq. Other kinds of regular- 
ization have been considered for similar purposes, for instance the 0— volume used in [21] or those 
used in [TniEl. 

If we knew that £ q was a monotone function, this result would yield immediately an equivalent 
for £ q (A) as A f oo by Karamata's Tauberian theorem. Unfortunately we don't know that it is 
a function neither that it is monotone (it is not the case in general). Nevertheless, we have the 
following result. 

Theorem 1.2. i) £ q is a continuous function on (0, oo). 

ii) Assume that G is of the form (1 i.ffl with g — g\. Assume furthermore that G is non trapping 
(see below) and that the principal symbols of Pq and P\ coincide outside the region {r > ro}, for 
some ro large enough, then we have the complete asymptotic expansion 

£ q (A)~A"/ 2 ^fo fe A- fc , AT+^o. (l.H) 

fc>0 

The coefficients bk can be deduced from ^l.lUQ and in particular &o = (27r) _,l w„ J x [di>o/ e ] q . 

Note that the condition on the principal symbol means that the terms of order 2 of Pi — Pq are 
supported near infinity. For instance P\ can be a perturbation of Pq with Pq = — <9 2 — e~ 2r A ff °(°) 
near infinity (the 'product case') or with Pq associated to a metric with constant curvature near 
infinity. The latter can be of special interest in view of the recent results of ■ 

We recall that G is a non trapping metric, if for any compact subset K of T*X \ there exists 
T K > such that ^{K) n K is empty for all |t| > T K , if $ is the geodesic flow on T*I. The 
results already obtained on R n let us hope that the non trapping condition in theorem ll.2l can be 
relaxed to get a Weyl formula, i.e. an equivalent for £ q . However this is an open question. 

We shall use methods of semi-classical analysis to prove these theorems and we will consider 

H t = h 2 P t = H a + eV, he (0, 1]. (1.12) 
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We will use the notations E e (fi) for the spectral projection of H t on (— oo, /i) and 

R e (z) = (H e -z)-\ US) = e- HH '' h . 

In theorem ll.il we will choose ft = t 1 / 2 and in theorem ll.2l we shall set ft = A -1 / 2 and consider 
the rescaled scattering phases £ q (/x, ft) associated to H , Hi, near the energy fj, = 1, since one has 
clearly £ q (^ft~ 2 ) = £ q Gu, ft). 

The non trapping condition will be used to show that 

/ ||(r)- M /(ffe)(7eW(r)- M || oo dt < Cmjmv (1-13) 
Jr 

for some M > 0, all / € Cq° supported close to 1 and uniformly w.r.t. ft 6 (0, fto] and e 6 [0, 1]. 
Such estimates are essentially well known under the non trapping condition. They have been 
proved in the Euclidean case by Robert-Tamura @E| and simplified by Gerard-Martinez [20] (see 
also HH1), using the theory of Mourre [SZj with a conjugate operator (see appendix [5J) defined as a 
suitable perturbation of generator of dilations r.hD r + hD r .r. As explained by Hislop and Froese 
in ^Hj , the generator of dilations does not fit the hyperbolic framework and they built explicitly 
another conjugate operator which makes Mourre theory applicable. In appendix [51 we sketch the 
proof leading to (|1 . 1 H|> by combination of the ideas of [201 ESI an d HBj > which is necessary since 
(11.13(1 is only known for fixed ft in the asymptotically hyperbolic case |18l H^j . 

Note finally that if l|1.13fl can be improved so that t p (r)~ M f (H e )U e (t) (r)~ M has a polynomial 
bound in ft -1 , then our method would show that £ q p ^ has a complete expansion, obtained by 
differentiating l|l.llfl . 

Acknowledgments: I want to thank Peter Hislop for helpful discussions as well as Gilles Carron 
and Didier Robert for their interest and useful remarks. 



2 The basic tools of the proof 

The purpose of this section, which is of pedagogic nature, is to describe the main tools of the proof 
of theorem II .21 The formulas that we are going to display hold for a much wider class of operators 
than those defined on asymptotically hyperbolic manifolds and we want to separate the general 
ideas, sketched in this section, from the specific analysis of the hyperbolic context given in sections 
El and H 



2.1 Representation formula of £ q 

Let us first assume that V, defined by (|1.12(l . is compactly supported so that the Birman-Krein 
spectral shift function £i^(fi,h) is well defined for the pair Ho,H e (this is essentially standard 
but is anyway a consequence of lemma I4.10fl . Then we can use the well known Birman-Solomyak 
formula [2) 



tv(f(H e )-f(H )) 



[ ti(f'(H se )eV) ds. (2.1) 
Jo 



The left hand side of 1(2. ljl is — (£i, e , /') so, with the notations of subsection 1 1.21 we get directly 

«'■•<"' h) - -iy^) d *=-I tt (lH ds ' <2 - 2) 



G 



since both sides of 1(2.2(1 have the same derivative (in the distributions sense w.r.t. fj,) and vanish 
for fi <C 0. For the general case, i.e. V non compactly supported, we consider the family 

V {K) = 9(nr)V9(Kr) : k g]0,1] (2.3) 

with 6 € 8 = 1 near 0, and thus we can define the associated family of spectral shift functions 
^[ K } . We obtain the following representation formula for £ q : 

Lemma 2.1. In the distributions sense on R M , we have 

f q (/i,A) = -Iim{tr(^-ywH (2.4) 



I \ djl 



q 



Here Ee is the spectral resolution associated with = Hq + eV [ K "> . 

The proof of this lemma (given in section |2J is not very hard and formally obvious from (|2.2|l 
and the definition of £ q . Formula H2.4(l leads obviously to the study of distributions of the form 

tr ^/(ffW)Fwj (2.5) 

with / £ Cq°, / = 1 close to 1. Note that such a function / can be added for free since we only 
consider \i close to 1. In order to simplify our notations, we drop the index k in the sequel but the 
reader must keep in mind that we work with a perturbation of the form ((2.3(1 . 

As usual, we consider the (semi-classical) Fourier transform of l|2.5|l which is tr (U e (t)f(H e )V). 
It is thus natural to consider distributions of the more general form 

tr (U e (t)f(H e )K e ) (2.6) 

with K e trace class, for instance K c = f(H e )V with / € C^°, // = /. Recall that we consider 

perturbations of the form (|2.3(l . so that f(H^)V^ is trace class for each k > 0, but for k = we 
only have K e £ S q in which case the trace makes sense only once {.} q has been taken. 

We will have to consider the e derivatives of f(H e )U e (t), that is why we quote quote the formula 

d e U e {t) = -J U e (t - s)VU £ (s) ds (2.7) 

which holds, for instance, in the strong sense on the domain of the operators H e (which is indepen- 
dent of e, see proposition ^. 9(1 . Since we want to use ((1.13(1 . it is important to keep a spectral cutoff 
in front of each U e (t) . To that end, we use a very simple trick. Let us introduce the notations 

U({t) = f(H e )U e (t), S f e = f(H e ). (2.8) 

Then for any small neighborhood I of supp /, we can choose / smooth, supported in /, such that 
ff — f and we obtain 

d e u{(t) = d e (s f e u e (t)s{) = (d e s{) u e (t)s{ + s{ (d e u £ {t)) s( + s f e u e {t)(d e s(). 

and using ((2.7(1 . we get the formula 

Mtf (.) - Mf um + i f um - .yvuh.) da + u'msl (2.9, 

This can be obviously iterated and proves the following more general result 
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Lemma 2.2. For all k > 1, d e U[(t) is a linear combination with universal coefficients of 

h~ j { d l °Sl Ul°{t )d\ x Si x V ■■■Vd l J- 1 S&- 1 Uh(t j )d l e iS{i d j (t ,- ■ ■ ,tj), (2.10) 

J Ft 

with l<j<k, l + --- + l j+ i = k-j, for-- JjJor-- Jj+i e Cg°(J) 
and of d e S(U( (t)d^ 1 S( , for <l < k. In ^2. KM . we have used the notations 

Fl = {(i , ■■■ ,tj)e [0, I to + • • • + tj = t} 

and dPifo, ■ ■ • ,tj) for the (j -dimensional) Lebesgue measure on the hyperplane to + • ■ • + tj =t. 

In the applications, we will get estimates on (|2.10|l using p. 13(1 combined with the following 
easy estimate 

dt < / |^o| * \i>i\ *•••* |^il (*) dt (2.11) 



Mto)Mti)---i>j(t 3 ) d 3 (t ,--- ,t 3 ) 

valid for all integrable functions tpor " ityj- 
2.2 Two microlocal tools 

The operators K c in (|2.6() will essentially be pseudo-differential operators. Sometimes we will need 
to shift the support of their symbols by the classical Hamilton flow <j)\. To that end, we will use 
the fact that for any to we have 

tr {U e {t)f{H e )K t ) = tr {U e (t)f(H e )U e (to)K e Ue(-t )) ■ (2.12) 

This remark was used by Robert in 4& and follows trivially by centrality of the trace. In general, 
we cannot obtain an explicit formula for U e (s)K e U e (—s) and rather get an approximation. This 
implies that we have to study an error term and this is why we display the following explicit 
formulas. If Kf is C 1 w.r.t. s, satisfying K® — K e , we have 

U e (t )K e U e (-t ) = Kl° + ^-JJ U e (to - s) (hi^K - [H e , K s e ]\ U e (s - to) ds (2.13) 

and this leads to the following exact formula 

tx{U e {t)f{H e )K e ) = tr (U e (t)f(H e )Ki°) + tr (u e (t)f(H e ) ° hid s K s e - [H e ,K s e ] ds^j (2.14) 

since the terms U € (to — s) and U e (s — to) cancel out by centrality. Thus, if we are able to find such 
a Kf with hid s Kl — [H e , Kf\ small in a certain sense, we see that the study of tr (U e (t)f(H e )K e ) 
reduces to the one of tr (U e (i)f(H c )Kl a ), up to a remainder which is given explicitly by (|2.14() . 
The method leading to the calculation of Kf is the usual one given by the Egorov theorem. We 
refer to for a proof of this theorem. The main point is that the symbol (in each chart) of K* 
has an explicit expression in term of the symbol of K e and of the Hamiltonian flow of p e . 

If K e is a pseudo-differential operator with a symbol supported in a suitable region of T*X 
(and this can be achieved by replacing K c by if* thanks to the above trick), we shall see that it 
can be factorized as 

K e = A e B* £1 with A e ,B, : L 2 (R n ) -> L 2 (X). (2.15) 
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In practice such a factorization will only be obtained approximately, i.e. K e — A e B* negligible 
(in a sense defined rigorously in section The operator A e will be used to intertwin U e (t) with 
a free dynamic given by U(t) — exp(—itP/h), i.e. make U e (t)A £ — AJJ(t) small, in some sense, 
with P = p(hD) differential operator with constant coefficients on R n . The explicit formula for 
U e (t)A e — A e U(t) is easily seen to be 

i r* 

U e (t)A e - A e U(t) = -rr / U e (t - s) {H t A e - A e P) U(s) ds (2.16) 

and this shows that, if 12. 1511 holds, then we have 

tr(U £ (t)f(H e )K £ ) -tr(f(H e )A e U(t)B*) = 

i r* 

- / tr (f(H e )U e {t - s) (H e A e - A e P) U{s)B*) ds. (2.17) 

Thus if we are able to show that (H e A e — A e P) U(s)B* is small (see lemma l5~^Ej) . we see that we are 
left with the study of tr (f(H e )A e U(t)B*). This trace is easy to study since we shall have explicit 
expressions for the operators A e ,B* and U{t). The construction of the operators A e and B e will 
follow the scheme of Isozaki-Kitada's method explained in the next subsection. 

2.3 The method of Isozaki-Kitada 

In this part, we recall the principle of the construction of the operators A e , B e by the method of 
Isozaki-Kitada introduced in [2S| • This method has only been used on M™ for Euclidean scattering 
[211 1461 1471 1151 121 but it turns out that it can be used in our framework as well, with some 
changes on which we shall put emphasize in section 

We first recall the algebraic formulas which enters into the game for a general differential 
operator of order 2 (on a general manifold X), which we still denote H t = h 2 P t (x,D x ) with 

P e (x, D x ) = Pe (x, D x ) + pM (x, D x ) + P W ( x ) 

in coordinates x — {x\, • • • , x n ). Here p e (x, £) is the principal symbol, i.e. homogeneous of degree 
2 w.r.t. to the dual coordinates £ = (£].,•■' ,£1), and p^ are homogeneous of degree 2 — j for 
j = 1, 2. We look for A e , B e defined as operators of the form J(ip e ,a e ) and J(ip e ,b e ) where 

J{<p t ,a e )u(x) = (2nh)- n J J e^^ x '^- x '^a e (x,^h)u{x') dx' df. (2.18) 

Note that l|2.18|l defines actually an operator from L 2 (M. n ) into itself (under suitable conditions on 
tp e and a e ); however, in the applications, a e (and b e ) will be supported into a region of K.™ x R^ 1 
whose projection onto R" is included into a coordinate chart of X. Thus, up to an invertible 
operator, we can consider (|2.18J1 as an operator from L 2 (R n ) to L 2 (X). 

Following i|2.16fl . we see that we have to study H e J(^p e7 a e ) — J(ip e , a e )P. Since P = p(hD) is a 
Fourier multiplier, we have obviously J(ip e , a e )P — J(ip e ,a e p). On the other hand, we see easily 
that 

H e (e*<*a E ) = (p e (x, d x (p e )a e + hi^L^x, d x )a e + h 2 P e (x, D x )a e ) . (2.19) 

Here, L e (x,d x ) is a differential operator of order 1 defined as 

L e (x,d x ) = w e (x,£).d x +c t {x,£) 
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where we have set 

w e (x,£) = (d^p e )(x,d x ip e ), c e (x,g) = p e (x,d x )cp e + ipP(x,d x ip e ). 
All this shows that if we look for a £ = ai°^ + ha^ + ■ ■ ■ + h N ai N \ then 

H e J(<p £ ,a e ) - J(<p e ,a e )P = J(<p e ,a e ) (2.20) 

with a e = ai 0) + ha^ H h h ( - N+2 ^a{ N+2) . For < j < N + 2, the functions a{ j) are given by 

5« = {p e {x,d x y e ) -p(0)aP -iL,{x,d x )a^- 1] + Pe(x, D x )a ( ^- 2 \ (2.21) 

(k) 

where we use the convention that a e — for k < or k > N. 

Since we want to make <|2.20[) small, we look for ip e (x, £) such that 

p e (x,d x <p e )=p($ (2.22) 

(k) 

which is usually called the Hamilton- Jacobi equation. We also need to find a\ such that 

L e (x,d x )ai 0) = 0, (2.23) 
L e (x,d x )a^ = -iP e (x ) D x )ai k - 1 \ k > 1, (2.24) 

which are the transport equations. The resolution of l|2.22|l . (|2.23|) and (|2.24l) rely upon estimates on 
classical trajectories. The technical part leading to such estimates in the asymptotically hyperbolic 
case is the purpose of section |3 Here we recall the general method. Assume that we can find 
S e (t, x, £) defined on [0, oo) x T, for some open set T C R 2 ™, such that 

d t S e =p e (x, d x S e ), S e (0,x,£) = x4- 

The existence of S e will follow from suitable estimates on <^>*, the Hamilton flow of p e . In practice, 
T is such that 

<9 5 S e (t,a:,£) -> oo, as tfoo. (2.25) 
Now, using the fact that S e is a generating function of the flow, i.e. 

4(x,d x S e ) = (d i S e ,0, (2.26) 
the invariance of p e by the flow and (|2.25|l imply that 

\imp e {x,d x S e )=p{0 (2.27) 

tjoo 

provided 

This shows that we have to built (p e such that d x tp e = lim^oo d x S e , or equivalently such that 

POO 

d xiPe (x,o = z+ / d t d x s e (t,x,o dt. 

Jo 

Of course, if such a function <p £ exists it is not unique. A possible construction is the following 

<p e (x,£) = x.^+ f°° d t (s e (t,x,0-Se(t,0) dt (2.28) 
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with S € (t,^) independent of x and such that the integral of the right hand side converges. In 
practice, S e will be easy to find. In that case, by construction, (|2.27|) holds and we get Q2.22[l , 

Note that, for Euclidean scattering on R™, we havep(£) = |£| 2 . In the asymptotically hyperbolic 
situation, with the coordinates x — (r, y) and £ = (p, rj) we will consider = p 2 . 

For the resolution of the transport equations l|2.23|l . (|2.24ll we study the solution x\{x, £) of 

dx l 

—±=w e {x\,£), x°{x,C)=x. 
Indeed, any solution a(x,£) of L e (x,d x )a = satisfies 

= a (x, £,) exp (- J c e (x s e , f) ds^j . 

In the applications, we shall prove that x\ is defined on [0, oo) x T and satisfies (in a sense to be 
made precise) 



x\ 



& + *Too. (2.29) 



This is actually well known in the Euclidean case (see for instance [291 1211 W!\ 115p . Hence, if we 
look for ac°' such that ai° (x, £) — > 1 as (x,£) — > oo in T, we obtain 

a^(x,0 = exp^ c e (xt,0dt). (2.30) 
Similarly, by the method of variation of constants, we see that any solution of L e a — i satisfies 

a{x\,0= ( a ( x '^ + ^ e^P ^ ^(^,0 duj ds^exp^-^ Ce ds 
Thus if we look for a solution at^ of (|2.24|1 such that ai (x, £) — * as (x, £) — > oo in T, we get 

ai k \x,0 = -j™iP e (x,D x )ai k -V(xt,0 exp c £ (i £ s , ds) df. (2.31) 

Here again, the convergence of the integrals is justified by the appropriate estimates on the functions 
c e and P e (x, D x )a{ k l ' which need to be shown. Such estimates are well known in the Euclidean 
case and will follow from section 03 for the hyperbolic one. 

The formulas l|2.30|l and l|2.31|) define ai k \ for k > 0, in T. In scction[5l we will explain how to 
define them globally, i.e. how to cut them off outside a suitable area. 

Recall that we want to consider a factorization of the form J(tp e , a € ) J(y c , 6 e )* whose Schwartz 
kernel is 

JC e {x,x') = {2Trh)- n [ e *((p.(*.0-»'.(« / ,0) Oe ( a ; ) £ h)b e {x', C, h) dC. 



This kernel is the one of a pseudo-differential operator since Kuranishi 's trick, namely 

¥>e[x,0 ~ <Pe(x',0 = (X- x').6 e (x,x' ,£) 

which is of course obtained by Taylor's formula, allows to write 

K, e (x,x') = (2nh)- n ( e^ {x - x ">- 9 d e (x,x',9,h) d9 
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provided the following map is a diffcomorphism for each x, x' belonging to the projections of the 
supports of a e and b e : 

t^9 e (x,x',0 = [ d x tp c (x' + t(x-x'),£) dt. 
Jo 

Note that, in view of (|2.28() . we see that if d x dtS e (t,x,£) is small (which will indeed be the case), 
then the map £ i— > 6 C is close to the identity and easily seen to be a diffeomorphism. In this case 
we have 

$ e {x,x',9,h) = a £ {x, 0~ l {x, x' , 9), h)b e (x', 97 1 {x, x', 9), h) |det d 6 9j l {x,x' ,Q)\ . 
Using a general elementary property of pseudo-differential operators, we have 

K e (x,x') ~ (2nh)- n Y^h j [ e^ x ~ x ' )£ ^ —ffiD%$ e {x, x' , 9, h) w=x ffl. 

3>0 \a\=j a ' 

By identification of the powers of h, this allows to find b e — & e ' + • • ■ + h N bi N ^ such that, modulo 
h , the right hand side is the expansion of the Schwartz kernel of K e , with the notation of (|2.15(1 . 
For instance, if the principal symbol of K e is cr(x, 6), we must have 

af\x i e- 1 {x 1 x i e))b { e\x 1 97 1 {x 1 x i e)) |det d e 9- l {x,x,6)\ = a(x,9) 
which implies that 

b{°\x, = a(x, 9 e {x, x, £))a (0) (x, O" 1 |det d^ e (x, x,Q\. (2.32) 

More generally, one can get explicit expressions for b^\ ■ ■ ■ ,bi and the important remark is 
that they are linear combinations of products of (derivatives of) (x, a e for k > 1, 

6 e {x,x',£) (evaluated at x' — x) and the symbols of K e evaluated at (x,9 e (x,x,£)). In practice, 
9 e (x,x,£) — £ will be small in the region that we will consider and T will be a neighborhood of 
supper, so that bi°\ 6c 1 ', • • • will only depend on a and on the values of a{°\ ai 1 ^, ■ ■ ■ and (p e in T. 
In particular l|2.30JI will ensure that ae ' doesn't vanish. 

3 Some estimates on the geodesies 

In this technical section, we prove long time estimates on classical trajectories, in suitable areas of 
T*X, needed to justify Isozaki-Kitada's method in the asymptotically hyperbolic case. 

3.1 Results of this section 

We consider the function defined for r £ K, p £ K and y,r) £ by 

Pe(r, y, p, rj) = P 2 + e~ 2r 9e(r, y, 77) 

where the function g e (r,y,r]) defines a metric on R™ -1 , i.e. g e is a smooth function which is an 
homogeneous polynomial of degree 2 w.r.t. rj and such that for some Co > 

Co'lvl 2 <9e(r,y,v)<C \v\ 2 , (3.1) 
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for all r G R, y, r\ G R n 1 and e G [0, 1]. We assume furthermore that g e is of the following form 

9e(r, y, rf) = g(y, rj) + ee- r g (e" r , y, rj) (3.2) 

where g and g satisfy the following estimates 

\%d2tfg{x,y,ri)\ + \d^g(y, V )\ < C^M 2 "^, * e [0,1], y, V G R"" 1 . (3.3) 

Note that in particular, we have <7 e (r, y, 77) — (/(y, 77) = 0(e~ r )\rj\ 2 . All these estimates are satisfied 
by the principal symbol of P e in any chart at infinity (R, 00) x O. Thus we can assume that this 
principal symbol is the restriction of some function p f satisfying the above estimates. We are 
going to prove estimates on trajectories of some vector fields and these trajectories will turn out 
to lye inside (R, 00) xflxl". This means that the results of the present section can be obviously 
considered as local results in T*X. 

Our first goal is the study of the Hamiltonian flow of the function p e , that is the solution of 

# = H e (4) 

where the notation ' stands for d/dt throughout the section and H e = H e (r, y, p, 77) is the Hamilto- 
nian vector field of p e , that is 



/ d p p e \ 


( 


2p \ 


d v Pe 




e 2r dr,g t {r, y,rj) 


-OrPe 


e" 


~ 2r (2ff e (r, y, rj) - d r g £ (r, y, 77)) 


\-9yPeJ 


V 


-e~ 2r d y ge{r,y,il) J 



(3.4) 



We denote the components of the flow by (r*, y\, p\, 77*). They are functions of the initial condition 
(r, y, p, 77). We remark that 4>\ is defined for t G R for if this was wrong, then \<f)\\ should blow up 
in finite time. This cannot happen since the conservation of energy, p e — p e o implies easily 
that <j>l is bounded, hence (j>\ — 4>° e = 0{t) can not blow up in finite time. 

Let us introduce some notations. We consider energy intervals of width 2w defined by 

/(if) = (1 — w, 1 + u>). 

In the end, w will be chosen small enough but for the time being we only assume that < w < 1/2. 
We also introduce the outgoing (resp. incoming) parameters 0+ > and <5 + > (resp. er_ > 
and 6- > 0) defined by 

a± = (1 =F w) 1 ^ 2 — S±. 

Of course, this makes sense provided < S± < (1 =F w) 1 ^ 2 . The last two parameters we shall need 
are a positive real number R > and an arbitrary open subset C R" _1 . We can now define the 
outgoing area T + (R, a +7 w, 0) and the incoming area T~(i?, ct_, w, Q) by 

T ± (i2,<T±,«; J fi) := {(r,y,p,v) e R 2 " | Po(r,y,p,v) € /(«;), ±p > -<J±, r > fl, y G tl} . (3.5) 

The point of considering such areas is that we have a splitting 

r+(R,o-+,w,n)ur-{R,o--,w,n) = {(r,y,p,ri) \ Po (r,y,p, V )eI(w) 7 r > R, y G 0} . (3.6) 

The first result is the following. 
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Proposition 3.1. Let < w < 1/2, a± > and f2 as above. There exists R large enough and 
C > 0, depending only on Co and a finite number of constants Ck. a ,0 in $3.3$ . such that the 
following estimates hold 

r* > r±t-C, (3.7) 
Pi > (3.8) 
\y\-y\ < Ce~ r (3.9) 

for all e G [0,1], (r,y,p,rj) G T ± (i?, cr±, w, fl) and ±t > 0. 

The meaning of these statements is that the properties are true on T + (resp. T~) for t > 
(resp. t < 0). This result says that the geodesies starting in outgoing (resp. incoming) areas stay 
in the neighborhood of infinity in X for t > (resp. t < 0). More precise estimates on the flow 
will be given in the theorem below. Note moreover that i|3. 9J) shows that y\ lye in any arbitrary 
small neighborhood of £1. It explains why our results can be localized in charts at infinity on X. 

Remark 1. This proposition, or rather its proof, shows in particular that the flow of p c in 
T ± (R, a±,w,fl) depends only on the values of p e in the region {r > R — C}. This implies that 
if we replace g by g^ depending on a parameter k G [0, 1], in H3.2J) . such that the constants Co 
and Ck l0 i,/3 can be chosen uniformly w.r.t. n for r large enough, then the above proposition is true 
uniformly w.r.t. k. In particular, this holds if one considers the principal symbol of l|2.3|l which 
involves 

~g^(e~ r ,y,r 1 ) = e(Kr) 2 ~g(e- r ,y,7 1 ). (3.10) 

Using this remark, we now claim that all the results listed in this subsection hold true uniformly 
w.r.t. k G [0, 1] if one considers H3.10(l . 

Theorem 3.2. Let w 1 a± and fl be as in vrovosition \3 . 1\ There exists R large enough such that 
for all 7 defined by = d^d l p d^d{, one can find C 7 > satisfying 



where p e = p e (r,y,p, 



<9 7 (r* - r - 2tp) 

|9 7 (yt-y) 
\d" (pI-p) 



< C 7 (t)e-V+V r \f,\ 

< C 7 e~ (j+1)r |??| 

< C 7 e- (j+1)r |??| 

< C 7 e-( j+1 >\rj\ 



(3.11) 
(3.12) 
(3.13) 
(3.14) 



for all e G [0, 1] provided the following condition holds 

±t>0 and (r,y,p,rj) G T ± (R,a±,w,n). (3.15) 

Corollary 3.3. For R large enough and for all (r,y,p,r)) G T ± (i?, a±, w, fi) the following limits 
exist for all e G [0, 1] 

yf = lim yl, r)f = lim rf e lim p\ = ±p] /2 

t — >±oc t—*±oo t—*±oo 

'/)• 
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This corollary follows very easily from the motion equations, proposition 13. II and theorem 13.21 
since, in particular, y* p\ and rf are OieT 21 ). 



Another important consequence of the theorem is the following. Since V Pi?) (p*, 77') is close to 
the identity matrix if e~ r r\ is small enough, we can hope that (p, 77) 1— » (p*, 77') is a diffeomorphism 
under suitable conditions. This is why we introduce 

V ± (R,e,w,S\) = {(r,y,p,rf) \ p (r,y,p,rj) 6 I(w), r > R, ±p > 0, !/£(], e- 2r g(y,rj) < e}. 

Before studying the above diffeomorphism, let us note that any area r ± can be reached in finite 
time from Y^. Precisely we have 

Lemma 3.4. Let < w < w < 1/2 and o~± > 0. Assume that SI is bounded and that SI <s SI' . 
Then, there exists R > large enough such that for any e > and R 1 > R, there exists T > such 
that 

4>{ (T ± {R,a ± ,W,Sl)) C ^(#,6,111', SI'), V ±i>T, ££[0,1]. 

This lemma follows easily from the previous results, the main tool being the fact that e~ r »rf e — > 
as ±t — > +00. Now let us consider the map 

$* : {r,y,p,rf) ^ (r,y,p*,»£). 

Then we have the following result. 

Proposition 3.5. Let Slo, Sl\ be bounded, connected and such that Sli <j= SIq. There exists Rq, Ri > 
large enough and Wo,Wi,Sq,£i > small enough such that for all ±t > 

r ± (J? 1 , e u Wl , Or) C (r+(.Ro, e , fio)) • (3-16) 

Furthermore, is a diffeomorphism from T (Ro, eojWOj^o) onto its range for all ±t > 0. //we 
denote the inverse map by (r,y,p,rf) t— ► (r, y, p*, 77*) we /lave 

|^(p* -p) I + -^j < C^e-W+^l^l, if d"i =d*d«d l p dPdi (3.17) 

wii/i C 7 independent of e € [0, 1], ±i > and (r, y, p, 77) € r ± (i?i, £1, iui, Oj ) . 

The main application of this proposition is the resolution of the following eikonal equation 

d t St=Pe(r,y,drSt,d v St), S?(0,r,y, p, 77) = rp + y.77 (3.18) 

with 5^ = Sf(t,r,y, p,rj,e) defined on r ± (i?i, si, w\, Sl{) for any ±t > 0. We can solve <|3.18() 
since any S which satisfy 

S(t,r,y,pl,ril,e) = rp + y.T] + tp e - {rd r p c + y-d y p c ) o cj> s e ds (3.19) 

Jo 

solves (|3.18|) . Thus the composition of the right hand side of (|3.19|l with the inverse of is a 
solution to H3.18fl . We shall show the following 

Proposition 3.6. We have a solution Sf 1 of ^S.lS^ on (R\ , s\ , W\ , Sli) for ±t > 0. Furthermore, 
if R\ is large enough and £\ small enough, we have 

\W (Sf -rp~ y.n T tp 2 ) \ < C 7 e-^ +1 >\ V \, dP = d^d l p d^d{ (3.20) 

for e e [0,1], ±t > and (r,y,p,rj) 6 r ± (i? i , e\, Wi, Sl x ). 
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Remark 2. Wc add this short remark in order to explain in what sense these constructions are 
continuous w.r.t. k, if we consider (|3.10l) . Let us recall again that all the coming proofs in the next 
subsections work uniformly w.r.t. the parameter k G [0, 1] since they rely on the estimates (13.111 
and ()3.3|l which are uniform w.r.t. k, when working with (|3.1Q(I . We will omit this parameter not 
to burden the notations but we hope that the proofs are explicit enough to make this uniformity 
clear. Then the continuity w.r.t. k G [0, 1], and actually the smoothness, is obvious since in the 
case of the flow (theorem I3.2[l it follows from the standard result of smoothness of O.D.E. with 
respect to parameters and in the case of the diffeomorphism ( propositions l3~Bl and l3~B)l it is due to 
the implicit functions theorem. 

We now give the proofs of these results. We shall only consider the outgoing situation since 
the incoming one can be treated similarly. 

3.2 Proof of proposition 13. ll and theorem 13.21 

Most of the results that we are going to prove rely on the fact that 

- d rPt = 2( Pe - p 2 ) + 0(e~ r )( Pe - p 2 ). (3.21) 
The first step is the following. 

Lemma 3.7. One can choose R large enough so that for all (r, y,p, rj) G Y + (i?, er + , w, tt) and all 
e G [0, 1] the following condition holds: 

r* 1 > R and p* 1 > for some t x > r\>R and p\> pi 1 , Vf > t x . 

Proof. First we remark that there exists R > and c > such that 

-d r Pe > cg(y,e- r V ), r > R, y, V £ M n_1 , e £ [0,1]. (3.22) 

This follows easily from l|3.1fl and 1|3.3|) . It is the first example of application of (|3.21|) . Then we 
consider the set 

T ={t>t 1 p s f > pi 1 and > R, Vs G [tx,t]} . 

It is clear that t\ G T thus T is a non empty interval. We shall prove that T := supT is +oo. We 
argue by contradiction and assume that T is finite. By continuity, it is clear that T G T, and that 
there exists T' > T such that p s e > pi 1 /2 for all s G [t\, T'}. Now using the fact r* = 2p\ we get 

K > rl 1 + (s - ti)p^ >R, h<s< T'. 

By (j*£23l, this implies that p% > on [t u T']. In particular p\ > p\ for s in a larger interval than 
[ii,T] which is a contradiction. □ 

In the next lemma, we show that reaches a positive value at some time t%. More precisely 
we show that p* 1 > 1/2 for some positive time t\ which is uniform with respect to the initial 
conditions in the outgoing area. 

Lemma 3.8. There exists R(w,5+) > and t\(w^S + ) > such that for all initial condition 
(r, y, p, rf) G T + (i?, er + , w, Q) and all e G [0, 1] 

r* 1 > R and p^ > i. 
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Proof. We first note that the result is clear if r > R and p > 1/2 by the previous lemma. Thus we 
assume that 5+ — (1 — w) 1//2 < p < 1/2. We also choose w < w < 1/2 such that w — w < 6+. By 
possibly increasing ^ we may assume that <|3.22[1 holds and that 

Po(r,y,p,rj) G I(w) and r > p e (r,y, p,ri) € I(w) V e e [0, 1] (3.23) 

since po — Pe — 0{e~ r ) on pjj" 1 (/(■!/;)). We may assume moreover that R is large enough so that 



Now we choose R such that 



< Y if r > and p e (r, y, p, rj) G J(to) for all e £ [0, 1]. (3.24) 



i?>i? + 4-^, with c+ = 2 1 / 2 - 1 - i 
c+(5 + 2 

We start by proving that for any (r, y, p, ?;) G T + (i?, cr+, w, f2) we have 

r* > i? and p' > p for all i>0, ££ [0, 1]. (3.25) 

To that end, we proceed as in the previous lemma. We consider the set 

T = {t > | r e s > R and p\ > p, V s £ [0, t}} 

which is non empty interval since G T and we show that T := supT is +oo. Assume that this 
is wrong, then T belongs to T. Note that we may assume that < on [0, T] since otherwise we 
can find t% G [0, T] such that p* 1 > and we get a contradiction using lemma l3~71 Then (|3.24|l and 
the conservation of energy yields 

pi > 2p\l\p\l 2 + p\) - 5 -± > c + 6 + (3.26) 

where the second inequality follows from the fact that p X J 2 > (1/2) 1 / 2 and 

p\ /2 + P > (1 - w) 1/2 - (1 - wf 2 +6 + >(l- 1/2^)5+. 
As a consequence we get p^ — p > Tc + S + and this implies easily that 



c+5+ 

since pf < and p > 5+ — (1 — w) 1 / 2 . On the other hand, the first equation of motion yields 

rj > r - 2TE+(w) 1/2 > R — 4T > R. 

Thus pT > p and > R. This shows that there exists T' > T such that p e s > po and r\ > R on 
[0, T') which is a contradiction and completes the proof of (|3.25|) . 

We shall now prove the existence of ti in the same spirit. Recall that we can assume that 
p < 1/2. Then there exists to depending on e and the initial conditions such that p\ a = 1/2 
otherwise p\ < 1/2 for all f > and 

pi > 2{p\' 2 + p\){p\/ 2 -Pt)- 6 f> 2(2^ - 1/2)(1 - 2-^ 2 )5 + - 5 -±> S f (3-27) 
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since we can use lj3~23|l thanks to (jXIZljl) . This implies that p\ > p + tS + /2 -> +00 as t — » 00 which 
is forbidden by the energy conservation. This is a contradiction thus p\ reaches the value 1/2 at 
some positive time. Moreover (|3.27|l holds as long as p\ is lower than 1/2, thus if to is the smallest 
positive time for which p* f ° = 1/2 we have 1/2 — p > toS + /2. This yields 

t < ^-(a+ + l/2). 

Using lemma (|3.7ll it is clear that pi > 1/2 for t > (2cr + + l)/5+ and this completes the proof. □ 

Proof of vrovosition UTli We choose R> R with R as in lemma Ij3.8|l . Then for t >ti, with t\ as 
in lemma I5~%1 we have 

r\ > (t-tO+rJ 1 . (3.28) 
On the other hand, for t G [0, ti] we have 

r* - r > -2cxt 

for any c\ > p e (r,y, p,rj) (note that such a c\ can be chosen uniformly with respect to e and the 
initial conditions in T + (R, a+, w, Q)). These two estimates yields l|H.7|l . The proof of H3.8|) follows 
directly from (|3.21(l since one knows that e~ r ' is small uniformly for t > 0. Finally, the conservation 
of energy implies that e _r «?7* is bounded and thus y\ = 0(e~ r ~ f ), by (|3.7|l . This implies l|3.9|l and 
completes the proof. □ 

The rest of this subsection is now devoted to the proof of theorem 13.21 We start with the 
following lemma 

Lemma 3.9. If R is large enough, there exists C > such that 

\vl-v\<Ce- r \ V \ 
for all (r, y, p, rj) G T + (R, a+, w, fi), e G [0, 1] and t > 0. 

Proof. Choosing R large enough and using proposition 13 . II shows the existence of c > such that 

\(dyg e ){tt,yi,4)\<cg e (tt,y\,Tit) 

since r* is large. The last motion equation and the conservation of energy then show that rf c is 
bounded and thus rf e — r\ = 0(t). Putting this new estimate into the last motion equation and 
using the fact that e _2r « = 0(e~ 2r ~ 2t ) show that 77* — 77 is bounded. Since e~ r rj is bounded on 
T + , the lemma is proved provided e~ r r\ is away from any neighborhood of 0. Thus we assume now 
that e~ r rj is small enough so that p > 0. Since pi > and p > we have 

p 2 <(pI) 2 <p, 

and then, using the fact that e~ 2r *g e (rl, yl, rf e ) — p e — (pi) 2 we obtain 

\r,l\<c(p e -p 2 ) <cCo{e- r H) 2 . 

This shows that 77* — r\ — 0(te^ r \rj\) and putting this estimate into the motion equation of 77* as 
before we obtain rf e — r\ = 0(eT T \r}\) thanks to the exponential decay in time of e _r ». □ 
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We explain the strategy of the proof of theorem 13 . 21 with the case d 1 = d e . By standard results 
on ordinary differential equations we know that <p\ is smooth with respect to e. In particular, if we 
consider the matrix 



( drpPe d^ p p t d^ p p, 
d r7 jP e dyriPe ^ptjP- 



-d; r p e -d* r p e 

32 



3 



dLpe \ 



-d^ y p e J 



(3.29) 



M e (r,y,p,r)) = 

\~@ryPe ^yyP e 

applying d e to the motion equations yields 

X\ = M £ (t)Xl + Y* 

with the notations X\ = 9,0*, M e (t) = M e (<^) and Y* = (<9 e H e ) {<t>f). By (gUJ and we have 

M e (t) = M + 0(e X p(-2r*)5 e (r* } 2/*,7j*)), (3.30) 
Yt = 0(exp(-3r*)« ?e (r*, 2/ *,7 ?£ t )) (3.31) 

where the matrix M is given by 

/0 2 0\ 

M= 0000 

0' 

\0 0/ 

In particular, it is very easy to check that M — and this implies that 

exp(tM) = 1 + tM, t e M. 
We do a change of unknown function by considering X\ — exp(— tM)X\. It satisfies the equation 

X\ = M e (t)Xt+Y* (3.32) 

where the matrix M e (t) and the vector Y e * are given by 

M e (t) = exp(-*M) (M e (t) - M) exp(tAf ) 
Y* = exp(-iM)r e t . 

Lemma 13.91 and proposition 13 . II imply easily the following key estimates 

M e {t) = 0((t) 2 e- 2t - r \ V \), (3.33) 
F e * = O ((t) exp(-2t - 2r)\r)\) . (3.34) 

We shall deduce the estimates on X\ from l|3.32|l . (|3.33(l . l|3.34|l and the well known GronwalPs 
lemma which we quote under the following form. 

Lemma 3.10 (Gronwall's lemma). Letu{t) > be a continuous function on [t.0,+00) for some 
to > 0. Assume that 

u(t) < a u(s) ds + b, t > to 
J t 

for some a > and b > 0. Then we have 

u(t) < be a{t ~ to \ t>t a . 
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Since X® = X® — 0, we can turn l|3.32[l into the following integral equation 

X* = f M e (s)X s e ds + f Y e s ds. 
Jo Jo 

Then GronwalPs lemma combined with the estimates (13.34(1 shows that for any T > 0, there exists 
Ct independent of the initial conditions in T + (R, a+, w, £1) and of e such that 

\Xl\ < C T e- 2r \ril < t < T. (3.35) 
We can also write (13.321) as follows 



Xl= [ At(s)X°ds+ [ Y e (s)ds + Xj, t>T. 
Jt Jo 

In particular, by choosing T large enough we can show easily that 

\Xl\< £\X^s)\d S + C^e- 2r \ V l t>T 

for some C' T > 0. This is a simple consequence of the fact that M e (t) — > as t — > +oo. Thus 
another application of Gronwall's lemma shows that for some C' T > 

I*' I < C^e- 2r Me (t - T) , t>T. (3.36) 
Putting l|3.35|l and l|3.36fl into l|3.32|l and using 13.33|) . one gets 

X* = O ((tfe-'e- 2 ^) , t>0. 

We can now come back to X\ since one checks easily that X\ = MX\ + exp(tM)X* which implies 



X\ = exp(iM) / X* ds. 
Jo 



Using the explicit form of exp(iM) we conclude that 

9V* = O ((*)e- 2r M) , (3.37) 
|0 6 ^| + |ay e | + \Frf e \ = O (e- 2r \ V \) . (3.38) 

This completes the proof of theorem 13. 21 if d' f — d e . The other estimates for I7I =0 can be proved 
easily using lemma l3~4l and the motion equations. For |7| 7^ 0, we proceed by induction by applying 
d 1 to the motion equations. □ 

Remark. Using the results of theorem 13. 21 we can improve the estimates on y\. Since 

Ve = V+ [ e ~ 2r! d v9e , yl , 77* ) ds 
Jo 

where e~ r ^ d v g e (r^ , y|, fy|) is bounded, we obtain easily 

-y)\< C^e-^ +1 >\r,\, j+ = max(j, !)• (3.39) 
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3.3 Proofs propositions [3~5l and 

We start with a technical lemma. 

Lemma 3.11. Assume that f2 is connected. Then T + (R,s,w,Q) is connected. 

Proof. Consider (tq, yo, po, rj ) and (r±, y±, pi, 771) in r + (_R, e, w, $1). We prove the result under the 
conditions rjo ^ and r\\ ^ (the other cases are simpler). We first connect (ro, yo, Po, ?7o) to 
(fij 2/1, Po, rr ?o) for some r > by considering 

((1 - s)r Q + sr 1 ,y(s),p ,T(s)ri Q ), s 6 [0, 1] 

with y(s) a path joining yo to yi in f2 and t(s) the unique positive continuous function such that 

exp (-2((1 - s)r + sr x )) g(y(s), t(s)ijq) = exp (-2r ) g(yo, ?7o)- 

Then we connect (ri, yi, po, t??o) to (ri, yi, pi, ?7i) using a path (p(s),r](s)) such that 

J7(0) = r?7o, 77(1) =771, Q(p(s),r](s)) G [£ ,#i], E j = Po(rj,Vj, pj,Vj), 5 = 1)2, 

where <2 is the quadratic form on 1" defined by Q(p, 77) = p 2 + e~ 2ri g(yi, 77). This is possible since 
the fact that n > 2 implies that the region defined by Q(p, 77) 6 [i?o, -Ei] and p > is connected in 
R™. □ 

We split the proof of proposition 13.51 into three lemmas. 

Lemma 3.12. Let < w < 1/2 and Q, C K™ -1 be an open subset. There exists e > small 
enough and R large enough such that is diffeomorphism from Y + (R, e, w, f2) onto its range, for 
all e e [0,1] and t > 0. 

Proof. The estimates Ij3.13l) and H3.14f) show that V$* is as close to the identity as we want by 
choosing e small enough. Under such a condition, is a local diffeomorphism onto its range, thus 
it is global iff it is injective. If <5*(r, y, p, 77) = $*(/, y', p', 77') then r = r' and y = y'. Moreover 
(|3.13|) and l|3.14[) implies that 

(p*, 77^) (r, y, p, 77) - (p*, 77*) (r, y, p', t?') = (p - p', V - 77') + 0(e) (3.40) 

thus if the left hand side vanishes we have p — p' = 0(e) and 77 — 7/ = 0(e). On the other hand, 
Taylor's formula to the second order combined with H3.13|l and (|3.14(l show that 

(p*, 77*) (r, y, p, 77) - (p*, 77*) (r, y, p', 7?') = (1 - 0(e))(p - p', 77 - 7/) + 0(e)|(p - p', 77 - t/)| 2 

This implies that 

I (pl,vt) (r,y,p,v) - (pl,vt) {r,y,p',v')\ = (1 - c(£))l(p - p',v-v')\ 

if the left hand side of Ij3.40|l vanishes and shows that p — p' and 77 = 77', if e is small enough. □ 

Lemma 3.13. Let Oi , O2 be connected open subsets o/R™ -1 such that fli <s ^2- Then for all £2 
small enough and R2 large enough there exists exists W\ 1 W2 1 £\ > small enough and i?i > large 
enough so that i\3.16} holds for all t > and all e <E [0, 1] . 



21 



Proof. The proof follows a rather standard scheme (see . It is enough to find parameters such 
that 

r + (R 1 ,s 1 ,w 1 ,Q 1 )n^ t e (dT + (R 2 ,e 2 ,w 2 ,n 2 )) = (3.41) 
T+{R ll e 1 ,w 1 ,n 1 )n^l(r+(R 2 ,e 2 ,w 2i n 2 )) ^0. (3.42) 

Actually (|3.42|l is always satisfied since any element of the form (r, y, 1,0), with r large, is invariant 
by <!>*. Thus we look for conditions on the parameters Rj, Sj, Wj, flj ensuring the fact that l|3~4*T) 
holds for all i > and e G [0, 1]. Any clement in the intersection of (|3.41|) can be written 

(r, V, P, V) = {r, V, pl(r, y, Po, %), vt(r, y, po, Vo)) 

with (r,y 7 p ,rio) £ dT + (R 2l e 2l w 2l £l 2 ). Since Qi (s £l 2) choosing Ri > R 2 implies that one of the 
following two conditions is satisfied 

P 2 o + e- 2r g(y,vo) = l±w 2 , (3.43) 

e- 2r 9(y,Vo) = e 2 . (3.44) 

Using l|3~T3|l and (|3~14) we see that 

p 2 +e- 2r g(y,r,) = p 2 + e"* g(y, m ) + O (e\' 2 ) 
= l±w 2 (l + o(e 2 /2 /w 2 )) • 
Thus (j3.43(l can certainly not happen if the following conditions are satisfied simultaneously 

wi < w 2 , e 2 < w 2 , (3.45) 
since p 2 + e~ 2r g(y, rj) = 1 + O (wi). Similarly, we deduce easily from theorem 13 . 21 that 

e' 2r g{y, rj) = e~ 2r g(y, Vo) + O (e^) = e 2 (l + {e- R -/e 2 )) . 
Thus the conditions e~ 2r g(y, r/) < e± and (|3.44(l cannot hold simultaneously if 

El < e 2 and e~ Rl < e 2 . (3.46) 
The existence of parameters satisfying (I3.45|) and l|3.46|l is clear. We choose can choose for instance 
w\ = w 2 , e 2 — w 2 E\ = u>2, e~ Rl — w 2 , R\ > R 2 . 

Thus (|3.41() and (|3.42() hold. Since T + (R i, £i, wi,fti) is connected and $* is an homeomorphism 
on a neighborhood of T + (R 2 , e 2 , w 2 , fl 2 ) onto its range (by possibly decreasing w 2 ,e 2 ,fl 2 and 
increasing R 2 ), this implies IpHfif . □ 

Lemma 3.14. If £\ > is small enough, the estimates \3.1 r /\ hold on T + (R±, £\, w±, tti) fort > 
and e € [0, 1]. 

Proof. By (|3.16|) we know that (r,y, p\,i)\) € T+(R 2 , e 2 , w 2l fl 2 ) if {r,y,p,rj) 6 r+(iii, e\, w u Qi). 
Moreover if W2 is small enough, then r + (i?2, £2, W2, f^) C T + (i?2, — 1/2, u>2, ^2)- Thus we can 
use the estimates (|3.13JI and Ij3.14|) with the initial conditions (r, y, p*(r, y, p, 77), ry*(r, y, p,rf)). This 
shows that 

|»7* -r?| <Ce- r \f,i\ 
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and by choosing r large enough so that Ce r < 1/2 we get 

\f) t e -r ] \<2Ca- r \r)\. 

Using this estimate and theorem !3.2l we obtain similarly 

\fc-p\<Ce- r \r,i\<Ce- r \r ] \. 

This shows that (f3~TT)l holds with 7 = 0. For |-y | > 1 we proceed by induction, by differentiating 
the equality $* o ($*) = id. For instance, if d 1 = d e , we have 

(v#*) l( ,j r i.0 e (i>*)" 1 = -(a e ^) K **)-i 

where the right hand side is 0(e~ r ) and, on the left hand side, (V<£>*) — 1 is small enough hence 

3 e (Qiy 1 = -(Vfct)- 1 .^**)" 1 - 0(e- r ). 
We don't go any further into details. □ 
This lemma completes the proof of proposition 13.51 

Proof of vrovosition \3.b\ The fact that S^~, defined as the composition of the right hand side of 
(|3.19|) with the inverse of <£•', solves (|3.18ll is a standard result. See for instance ^5] or 03- Thus 
we focus on the proof of H3.2()(l . We first remark that H3.19[l can be rewritten as follows 

r\p\ + yUl - 2 f (p s e ) 2 ds + t Pe - [ e-<d n g e {rl,yt, if t ).e~<rf e ds. (3.47) 
Jo Jo 

This is easily obtained by integration by part using the motion equations. Furthermore 

rt 



f {plf ds =tp e - f e- 2r 'g t (r s e ,y s e ,v s e) ds 
Jo Jo 



and this implies that 

EH = r\p\ + vl.rt - tp e - [ e- r °d n g e (r s e ,y*, nD-e^H ~ 2e~ 2 < g^yl^ds. (3.48) 

Jo 

where one must notice that the integral is convergent and (D(e~ r \r)\). We rewrite ({3.48(1 using the 
fact that p e = p e o 0*. In particular, if we note ff = r|(r, y, p', ff t ) and y s c — y|(r, y, p\, 77*) we obtain 

S+{t,r,y,p,r],e) = v\p + y\.r] - t [p 2 + e- 2f 'g e {rl y\, 77)) + O (e^l) 

= flp + y.ri-tp 2 + 0(e- r \r]\) 

using l|3.39|l to estimate y\ — y and the exponential decay w.r.t. t of e -2r «g e (f*, y*, 77). In order to 
estimate f*, we use the motion equations to write 

rl = r + 2 f pl{r,y,p\,fi\) ds, p s e (r, y, p\, fjt) = p + f (d rPe ) (f-^y^^i^y, pl^)) du. 



By (|3.7|) and theorem l3.2l we see that the second integral is 0(e s e r \?]\) and this implies that 

f\ = r + 2tp + 0{e- r \r 1 \). 

and finally we obtain 

S?(t, r, y, p, r], e) = rp + y.r) + tp 2 + 0(e~ r |?7|). 

Note that the remainder 0(e~ r \r)\) is explicit and theorem 13.21 combined with proposition ^ . 51 show 
that (|3.20() holds. This completes the proof of proposition 13. 61 □ 
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4 Pseudo-differential operators 
4.1 Local theory 

In this part, we consider pseudo-differential operators on R n . We apply the results to operators 
on X in the next subsection. 

Definition 4.1. For any m,m' £ R, the class S m ' m is the set of smooth functions a such that, 
for all R > 0, Q <s R™ -1 and k, I, a, (3 there exists C such that 

\d k r d l p d^d^a(r,y,p,r,)\ < C {p) m {e^ rj) m ' , 
for allr>R, yen, p £ R, n £ R n_1 . FFe set 5^°° = n m , m -5 m - m '. 

As usual, the best constants C are semi-norms which define the topology of S m ' m . We also 
mention that we shall mainly consider cases where m,m! € R~. 

We give two examples of special interest for us. If h belongs to S™{ m , that is if 

\d k r d l p d^b{r,y,p,ri)\ < C{p) m {n) m '-^ 

for r > 0, y e ft, pGl and 77 £ R n_1 , then one checks easily that the function a defined by 

a(r,y,P,v) = b(r,y,p,e~ r n) (4.1) 

is an element of S m ' m . In particular, if / £ S then f {p 2 + e~ 2r g(y, 77)) £ S* -00 . The second 
example is the following. If £lo d ft and i// > u>, there exists C > such that for all R large 
enough and all a £ suppportcd in T + (R, <r + , w, fl ) (see l|3.5ll ') we have 

a o cj)- 1 £ S^ 00 and supp a o $7* c T+(i? + i - C, o+, it;', O), (4.2) 

for all t > and e £ [0, 1]. This follows from proposition 13 . II and theorem 13. 21 
The pseudo-differential operators that we will use are of the form Oph(a), with 

Op h (a)u(r, y) = (2tt)-" / / e ir P +i ^a(r, y, hp, hn)u(p, 77) dpdn, h £ (0, 1], (4.3) 



where u(p, 77) — JJ e~ jrp ~ 4a ' ?, ii(r, p) drdy is the Fourier transform of u £ S . They depend on the 
parameter h £ (0, 1] and it is natural to consider symbols depending on h as well. Following the 
standard definitions of 05], we say that a — a(h) is an admissible symbol in S m ' m and note 

a ~ do + hai + h 2 a2 + • • • . 

to mean that for all N, a = ao + ha± + ■ ■ ■ + /i JV_1 ajv_i + h N r^(h) with a,j £ S m,m , independent 
of h, and r N (h) bounded family of S m - m . 

We now give estimates on Oph(a) in Schatten classes. 

Proposition 4.2. Let \be a bounded function supported in R + x il, with f2 <e R™ _1 and v > 0, 
q > 1 be positive real numbers such that v > (n— l)/q. Then for any a £ S~ £ ~ 1 / q ~ e ~ tyn ~ 1 ^ q , with 
e > 0, the operator e~ ur xOph(a) belongs to S q and 

||e-- X ^(a)|| q <C/ l -"/ 9 , h £ (0,1], (4.4) 

where C depends on finitely many semi-norms of a. In particular, if ' v > (n— 1), then e~ ur xOph{a) 
is trace class and 

tr (e-» r xOp h (a)) = (2nh)- n [[[[ e — x (r, y)a(r, y, p, V )dpd V drdy . (4.5) 
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Proof. With no loss of generality we can assume that e is a small as we want and in particular that 

Ti — 1 

v>2e^ . (4.6) 

q 

We can also assume that x(r, y) = cf>(r)ip(y) with if) e C^°(fl) and cj> = l near infinity. The estimate 
(14.4|> will follow from the fact that we can write e~ vr xOph{a) = AhBh with 

\\A h \\^<Gh- n l\ \\B h \U <C, he (0,1]- (4.7) 

In order to construct Ah, we pick ip £ Cq° such that tpip = ip, and we set Ah — A\ ® A\ with 

A\ = e —0(r)(^ r .)- e - 1 /«, A \ = %{y){hD v )-°-W*. 

By standard estimates, we know that A\ = 0(h~ 1/q ) in S q (L 2 (tt, dr)) and that A 2 , = 0(/i-("- l )/«) 
in Sq(L 2 (R n ~ 1 ,dy)), thus the first estimate of l|4.7|) holds. We now have to consider Bh which we 
define, using <fi such that (jxj) = </>, by 

B h = {hD r y +1 l^hD y y+( n ^e-^ r 4>{r)^{y)Oph{a). 

In order to show the second estimate of H4.7(l . it is enough to show that Bh = Oph(b) with b = b(h) 
bounded in 5°'°. We first consider {hD y ) e+ ^ n -^ q e-^- e >^{r)i){y)Op h (a) = Op h {bi) with 

6i(r, y, V, /i) = e-^^ir^ir) 1 -'" J e^(hr} + h^f+^'^r, £, ftp, ^)cie (4.8) 

where a,/, is the Fourier transform of ip(y)a(r, y, p, ?y) with respect to y. By Peetre's inequality we 
have (hrj + ft^+C™- 1 )/? < CXftr/}^™- 1 )/^^}^™- 1 )^ and the integrand of gSJ is dominated, 
for any M , by 

C(e) _M e-^- e )'"(/i7 ? ) e+ ( n - :L 5/5(ftp)- e - ;L /«(/ie- r 7 ? )- e -( n - 1 )/« < C(0 _M (V) _£ " 1/<z , 

where we used (14.611 . The same holds for the derivatives of b\ and shows this b\ 6 S* -6-1 / 9,0 . 
Furthermore it depends continuously on b. Thus (hD r ) E+1 / q Oph(bi) — Oph(b) with & bounded and 
the second estimate of (14 .7|) holds. The proof of (|4.5(l is well known. □ 

This proposition is the key of theorem ll.il . For theorem 1 1.21 we shall need the slightly stronger 
estimate H4.1()[) below, in order to use H1.13|) . By standard results on pseudo-differential calculus 
we know that, for any MsR and b € S°'° 

(r) M Op h (b)(r)- M = Ophib), b e S°'°. (4.9) 

Furthermore, for all M, (r) M e~ er <p{r){hD r )- e - 1 /i (r) M still belongs to S q with norm 0(ft~ n /«). 
Using this remark and the proof of proposition 14. 21 we see easily that 

\\ e -^(r} M xOp h (a){r} M \\ <Ch- n/q , h e (0,1], (4.10) 

for any M £ R, a 6 S~ e ~ 1 / q '~ £ ~( n ~ 1 ^ q and v > (n — l)/q. Furthermore, if a (K) is a family of 
symbols which is bounded in £-e-i/<r.-e-(n-i)/g such that a («) a (o) in c°o( R 2n) ( or p'^n)^ 

as k J, and such that p 2 + e^ 2r \rj\ 2 is bounded, independently of k, on their support then 

^ w (r) M X*(« W )W M -^"« M X*(« ( ")W M mS q , « 1 0. (4.11) 
This follows simply from the fact that e^ v a (K) -> e~ d> a (0) in ^-e-i/g ,-e-(n-i)/g for any s > Qj 
and from the fact that we can choose 5 small enough such that v — 5 > (n — l)/q. 

Now we are going to study symbols and operator depending smoothly on e £ [0, 1]. 
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Definition 4.3. The class S" ,m ' m is the set of functions a e (r,y, p,rj) which are smooth w.r.t. 
r, y, p, rj and e such that for all R > 0, 51 <s R" _1 and j, k, I, a, f3 there exists C satisfying 

\did^d l p d^a e (r,y,p,r,)\ < Ce~^ +V > (p) m -(e^r,)™' ', 

for allr>R, y £ CI, pel, rj £ R™ -1 . We set S^'°° = C\ m , m ' S»> m > m ' . 

The typical example of such symbols is given by a e (r,y, p,rj) = f(p 2 + e~ 2r g e (r, y, 77)), which 
belongs to a e £ S°'~°° if / € S. If we replace a by a e £ S*, -°° then (JOJ holds with 5° ~°° instead 
of S~°°. This class is particularly natural and convenient since proposition 14. 21 shows that for all 

j \d'i xOph{ae)\ | q/j < Ch' nj/q , provided ?>n-l, g>j>l. (4.12) 

The main drawback of these classes is the following. If a € £ gv,m,rn ^ nen we don't have in 
general Qj/j(a e )* = Oph(a e ) with a c £ S^ m > m , which is due to the fact that pseudo-differential 
operators do not preserve exponential decay. This problem can be overcome by considering properly 
supported operators. This is the purpose of what follows. 

Recall that the Schwartz kernel of Oph{a c ) is given by the following oscillatory integral 

)C h , e (r,r',y,y') = (2nh)- n JJ S^^^^a^y, P ,rf] dp d v 

which we can write, using 9 £ Cq°(R) and 9 = 1 on (—6, 5), as K,h, t — K-hV + ^he with 

/C*f(r, r', y, y') = K h , e (r, r', y, y')9{r - r')9{\y -y'\). (4.13) 

Definition 4.4. If a e £ S" ,m ' m , we define Op p h r (a e ) as the operator with Schwartz kernel ICf™ 9 . 

The following proposition explains in which sense Op^ (a e ) — Oph{a t ) is small. 

Proposition 4.5. Assume that v > 0, m < and ml < 0. Let x be a bounded function supported 
in R + x ft, with fi <s R n_1 . Let Kh.t the operator with kernel x{ T i y)^"h e( r ' r '> V> f')* Then for all 
s £ R, j > and N > 0, there exists C > such that 

<Ch N , h £ (0,1], e£[0,l]. (4.14) 



^ + » )r (r)- s —K Kt (rY 



Proof. It is standard. On the support of /C^f{, we have either r — r'\ > 1 or \y — y'\ > 1 and we 
can do as many integrations by parts as we want with \r' — r\ hd p or \y' — y\~ 2 h 2 A ?; . This shows 

-off I 



that d^lC? (r,r',y,y') is a linear combination of 

h N (2nh)- n e- jr { { e i ( r - r '>/ h+i ( y - y ')- r >/ h b N , e (r, r' , y, y', hp, hrj) dpdrj (4.15) 



with &jv,e a product of derivatives of d{a e , x an d \r — r'\ 1 , y — The result follows from (|A.1|I 
and pO|l . □ 

Note that the adjoint e satisfies (|4.14|) with e^ +,y) '' on the right. This leads to the following 
definition. 
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Definition 4.6. A family of bounded operators A €; h is said h N negligible if for all j > and 

M e R, we can write d 3 e A eyh = A° e h + A) h H h A j e h where, for < k < j, 

\\ e ^r {r) -M A k he u k r {r) M^ < C NM , 3 h N , h G (0,1], 6 £ [0,1], 

for some pair fk, ^j-k > such that i/j + Vj-k = v + j ■ It is called negligible if it is h negligible 
for all N. 

Proposition 14.5(1 shows that xQPh ( a e) — xQphifle) is negligible, provided m < and m! < 0. If 
A c j l is negligible, then A* h is clearly negligible as well, hence all this shows that, by using standard 
methods for the calculus of the adjoint of a pseudo-differential, we obtain easily the following result. 

Proposition 4.7. Let a e £ S^ m,m be supported in (0, oo) x Q, with Q bounded and m,m' < 0. 
Then Oph{a e )* is the sum of Oph(a e ) and of a negligible operator, with a e S gv,m,m suc /j 

We end this subsection with a simple remark. In practice, we will use the negligibility as follows: 
if A c ,h is h negligible, we can write for all j > 1, M > and e > as 

j 

dtA e , h =h N J2 e~" kr {r)- M B* h (r)- M e^- k (4.16) 
k=a 

where the operators h are bounded in operator norm and Vk + Vj-k = v-\-j — e, with Vk, Vj-k > 
for all k. 

4.2 Global theory and functional calculus 

The pseudo-differential operators that we are going to consider on X will be of the form 

A = £ XkOph(a k ) X k, I finite (4.17) 
kei 

For each k, \k and \k are supported in the same chart which is either relatively compact or a chart 
at infinity, and a k is a symbol expressed in the coordinates associated to the chart. If we work in 
a chart at infinity, we will always use the radial variable r and variables y±, ■ ■ ■ , y n -\ associated to 
the manifold at infinity Y . 

Let us recall some (standard) abuse of notations which are convenient. We use the same notation 
for XkOph(ak)xk as an operator acting on L 2 (W l ) or on L 2 (X). Furthermore if, for each k S J, 
XkOph(a k )Xk is bounded on L 2 (M™) with the Lebesgue measure associated to the corresponding 
coordinates, then A is bounded on L 2 (X) and its norm can be estimated by the sum of norms of 
XkOp h {a k ) X k on L 2 (R"). This is due to our choice of the density dvol. The same remark holds for 
estimates in Sq(L 2 (X)) and we will therefore use the notations ||.||oo and ||.|| q , initially used on 
i 2 (R"), for the respective norms of bounded operators and Schatten classes relative to L 2 (X). 

We now apply the results of the previous subsection to the analysis of functions of H t . We will 
only consider Schwartz functions / and use the Hclffer-Sjostrand formula, 

f(H e ) = -!- // dj(s + it)R e (s + it) ds dt (4.18) 

^ J JR2 
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where dz = d s + idt and a / is a quasi-analytic extension of /. i.e. a C°° function on C, supported 
in the strip |t| < 1 and such that 



/|R = /, 



M I \-M 



d s f(a + it) <C M \t\ M (s) 



for all s, t G R and M > 0. We do not insist on the construction of / which can be chosen 
depending continuously on / and refer for instance to ^B] for the details. 

The formula (|4.18(l shows that we only have to study the pseudo-differential expansion of R e {z). 
We recall that the pseudo-differential analysis of the resolvent is well known for operators on 
compact manifolds (or for elliptic operators on K") and thus we will only focus on the calculations 
in charts at infinity. We look for a parametrix Q e (z) of R e (z) of the form (|4.17l) and more precisely 



kei 



where J2 keI Xk — 1 is an admissible partition of unit, Xk = 1 near the support of Xk, and q* (z) is 
an admissible symbol for each k. We explain the construction of these symbols in a single chart at 
infinity and drop the index k for convenience. Since we want to get (H e — z)Q e (z) — 1 + 0(h°°), 
we seek q t {z) ~ . h? qj^(z) satisfying, in the chart that we consider, 



j,l,m 



h i+l+m (pf #3 . eW ) „ 1. (4.19) 



Here pf are the symbols of H e — z, i.e. Pq e — p e — z and pi tlL — pi for I = 1, 2, and the notation 
(a=ff=b) m stands for the m-th symbol of the (finite) expansion of the product Oph{a)Oph{b) if Oph{a) 
is a differential operator. The condition (|4.19|l yields 

q , e (z) = (p 2 +g e (r,y,e~ r r 1 )~z)- 1 (4.20) 

qj , e (z) = -qo, e V (p¥ o) #Qii,e) . J>1- ( 4 - 21 ) 



ja+ji+32=j, ji<j 



This procedure is the standard one used by Seeley [SOj and Helffer-Robert [22] but the point that 
we want to make here is the following: since the symbols of H e are of the form 1)4. it follows 
clearly that qj,e(z) defined as above is of the form l|4.1|) as well. By an easy induction, we get, for 
3 > 1, 

9i,e=^2dt Jte {p e -z)- 1 - 1 , (4.22) 
i 

where d^k,t is polynomial w.r.t. the variables p and e~ r rj which is independent of z and a linear 
combination of elements of S® ,m ' m with m + m' < 21 —j. In particular we can write 5j, e (r, y, p, Tj) = 
qj^ e (r,y, p, e~ r rf) for some qj >t € S Q \ 1 ^ 2 and we have 

q jie 6 S°' m ' m \ for all m, m' < such that m + m' = -2-j,. 

These remarks, combined with the fact that \p(p — z)~ x \ < C(Re z)/\Im z\, for p e R and z £ R 
show that the semi-norms of <7j ie in S® ,m ' m (m + to' = — 2 — j) are dominated by 



, (Rc z) M 
llm z\ 
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for some C, M depending on the semi-norm. 

We can now construct the global parametrix. We fix N > and define it as 

Q N , e (z) = ]Txfc I J2 h J Op h (qU z )) I Xk, 
kei \j<N J 

with the q\J e (z) defined in each chart by the preceeding procedure. Then we have 
(H e - z)Q N , e (z) = *M - z) E hj °Ph(qU z )) I U 



+ J2 [He,Xk] E hiOp h (ql(z)) ) Xk . (4.24) 

kel \j<N 

Lemma 4.8. For each k and j, the Schwartz kernel of [H e ,Xk]Oph(qj e (z))xk is 0(h°°) in the 
Schwartz space S(M. 2n ). More precisely, if we note it /C(r, r' , y, y' , e, z, h) then we have 

d? M Mry,y,y',e,z,h)\ < ch M e- Mr (r')- M (y)- M (y')- M ^^ 

for all multiindex j, M > 0, r, r' > 0, y, y' £ R" _1 , e G [0, 1], h G [0, 1] and z f R. 

Proof. It follows easily by integrations by parts similar to those of lemma 14.51 The exponential 
decay is due to the following fact. We can choose Xk = </'fc( r )V'/c(2/)j with <pk = 1 near infinity. 
Then [H £ , X k] is either compactly supported w.r.t. r, or \y — y'\ ^ in which case we integrate by 
part using A,, and get as many powers of e~ r as we want. We omit the other details. □ 

This lemma shows in particular that that the second term of the right hand side of H4.24f is 
negligible, in the sense of definition 14.61 On the other hand, by construction, the first sum of the 
right hand side of l|4.24|l is 

+ h N Op h (g h N ^z))) Xk = l + h N Y J XkOp h (g k N ^z)) X k 

kel fcGJ 

where Q% e (z) = Ylj+i+ m >N ^ + m (Pf,e#^i,e)m- It is not hard to check that it belongs to 



Se' N ^ 2 (for instance) using the form of q k e . We can summarize our result as follows. 



Proposition 4.9. The operators H e are essentially self-adjoint from C^°(X) and the domain of 
their self-adjoint realizations is independent on e £ [0, 1] . We have for all N 

Re(z) = QnAz) + h N R e (z)K N , e (z). (4.25) 

Here h TZ,N )fL (z) is the sum of^2 keI h N XkOph(g% e (z))xk and of the operators with kernel studied 
in lemma \Jl^ 

Proof. By the standard trick (see @H]), we see that (H e ± i)* is injective for h small enough since 
(H e ± «)<3iv,e(±i) = 1 + 0{h N ) in operator norm on L 2 (X). This implies the existence of a unique 
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self-adjoint realization. Hence the resolvent R e (z) is well defined and we can apply it to the left of 
(|4.24() which yields 14.251) . In particular, for h small enough and fixed z we have 

Re(z) = QnAz) (1 - h N TZ N , e {z))~ 1 

and since H tl Qiq t(i2 {z) is bounded for all £1,62, the same holds for H ei R e2 (z). This proves the 
independence of the domain w.r.t. e. □ 

Let us give a first application of this proposition. 

Lemma 4.10. Let A be a differential operator on X of order m of the following form in any chart 
at infinity 

A= a Ua (r,y)(e- r D y ) a D l r 

l-\-\a\<m 

with ai a bounded. Then for all is, q, k such that v > (n ~ l)/q and 2k — m > n/q , there exists C 
such that 

\\e- vr h m AR t (i) 1+k \\^<Ch- n ' q , he (0,1], ee[0,l]. 
If v = 0, by convention the above norm is the operator norm ||.||oo and h n l q = 1. 
Note that this lemma implies the following estimate, for Im z ^ 0, 

We-^AReiz)^]^ < ChT n '« r; l +k , he (0,1], ee [0,1]. 

Proof. We use the fact that k\R e (z) k+1 = d k R e (z) and Leibnitz rule into (|4.25|l to obtain 



RS) k+1 = ^ {d k z Q N ,e) (») ^1 - h N £ (H e - i)M^ |e ( < ) _j __ j 



(4.26) 



for h < ho small enough. Note that (iJ e — %Y dglZjy ,e{i) is bounded for all fi. Using the form 
of the symbols of Qiv.e given by (|4.22|) and by proposition 14 . 21 applied to e~" r h m Ad k QN,e{i), we 
obtain the result for h < ho. Note that we use the fact that k\d k qo^ c (z) — qo^(z) k+1 belongs to 

Se ' J £ ' 6 for all j, |a| such that j + \a\ = m and e > small enough. The result 

for h < 1 follows by writing i? e (i) = (/j 2 P e + i) _1 (/i 2 P e + i)R t (i) with h = mm(h , h). □ 

This proposition will be used in order to estimate d{(R,:(z)TlN,e(z)) in Schatten classes. We 
shall also need estimates similar to (|4.10() . To that end, we first recall that for any < Mo < 1 

[(r) M °,R e (z)] = -R £ (z) [(r) M °,H e ] R £ (z) 

where [(r) M °,i7 e ] is a differential operator of order 1 which is H e bounded since [d r ,(r) Ma ] is 
bounded. This implies in particular that 

{r) M «R e (z){r)- M <> = R c (z)~R e (z) [<r) M ° , H e ] R e (z){r)- M ". 

Since any M > can be written IMo with I £ N, this formula can be iterated to show that 

,(Re z) M ' 



\{H t + i)(r) M R t {z)(r)- M \\ <C- 



Im z 



M' 
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for some M' > 0, uniformly w.r.t. h £ (0, 1] and e £ [0, 1]. More generally, we can obtain rather 
easily for any k £ N and M > 

||(fl. + 1+fc (r>"ii e W 1 + fc (r>- 4 '|| oo <C^^. (4.27) 

for some C and M' independent of e and ft,. On the other hand, d{R t {z) = (— i? e (z)(V"i? e (z)) : ' 
can be written for any MeMas 



^i? e (z) = (-l)'j! (r)- M «r> w i2e(«)(r>- M ) (F(r)- M i2 e (^)(r) M ) J (r)" 



where the differential operator V — (r) M V(r) M is of the form e vr h 2 A for any v < 1 with the 
notations of lemma El. 101 Using these remarks combined with lemma H. 101 we will obtain 

Lemma 4.11. For any q > n — 1, M > real and j > 1 integer, there exists C, M' such that 

\\(r) M (diRM) W-Vril^ < Ch-^^^r 
for all h £ (0, 1] anc! e £ [0, 1] . 

Proof. We proceed by induction on j. If j = 1, we write 

(r) M d ( R e (z)(r) M = ((r) M R e (z)(r)- M e-» r ) (h 2 A( r y M R e (z)(r) M ) 

with (n — l)/q < v < 1 and e~ ur h 2 A = (r) V(r) . The first factor belongs to S q whereas the 
second one is bounded, by lemma l4~THl If j > 2, we have d(R € (z) = —j(d^~ 1 R e (z))VR e (z) and 
then the operator that we want to estimate can be written 

-3 {(r) M {dl^Reizj) R e {i) j - 2 (r) M ) {(r)- M \H t + iy- 2 VR £ (z)R e (i)(r) M ) • 

We use the the induction assumption for the first factor, which belongs to Sq/Q_i). We estimate 
the second factor using lemma l4~TOl and (|4.27|) and conclude using (jl.5|) . □ 

The main consequence of this lemma is the following. 

Proposition 4.12. For each N > No large enough and M > 0, q > n — 1 as above and j > 1 
there exists C and M' such that 

\\{r) M di (R e (z) - Q N , e (z)) (r) M || q/ . < Ch N - n ^ { ^™ (4.28) 

for all e € [0, 1] and h £ (0, 1]. 

Proof. Since h N d{ (R e (z)lZjsi ,e(z)) which is a linear combination of h N di 1 R e (z)dl 2 lZ]sr t e(z) with 
ji + ]2 — j, we have to estimates operators of the form 

h N ^ r )M {d n Rt{z))Re{l) n-i {r) M^ ( (r )-«( ff£ + i )ii-l^^ £ ( z )( r )M) _ 

The first factor can be estimated by proposition 14. 1 II and the second one by l|4.10|l and lemma 
□ 

Using formula (|4.18() . we get directly the following result. 
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Theorem 4.13. Let q > n — 1, / G S(M) and N > 1. There exists a pseudo-differential operator 
Q N e = X)/<at ^'^( iflii/i Ai of the form ( |^._?7| ) ; mf/i symbols in S® swc/i i/iai /or a/Z A/ > 

< C*^^'/ 9 , ft G (0,1], e G [0,1]. 

q/j 

_ffere C depends on a finite number of semi-norms of f. In each chart, the symbols of Ai are linear 
combinations of di'j lC f^ (p £ ) and in particular the principal symbol is /(p c ). 

4.3 Proofs of theorem 11.11 and lemma f2.ll 

Theorem 11.11 is a direct consequence of theorem 14.131 We obtain (|1.1U|) by considering h = t x l 2 
and / G iS(IR) such that /(A) = e~ A near the spectra of the operators. A priori, (|1.10|) involves all 
powers of the form £( fe ~™)/ 2 7 but a standard argument shows that the coefficients corresponding to 
odd k vanish, since they correspond to integrals of odd functions on the sphere. □ 

For the proof of lemma IP we have to show that tr{/(-ffi K) )FW} q -> ty{f(H e )V}q as k | 0, 
for all / G iS(K). Using the explicit expressions of the symbols of associated to (with 

the notation of theorem l4.13|) . it is easy to check that 

{v (k) QnY^ -> [vQ f Ne j in the trace class 

with a trace norm uniformly bounded by C(Re z) M /|Im z\ M for some M and C. Thus we are left 
to the study of the remainder which, via Helffer-Sjostrand formula, reduces to the study of the 
remainder given by (|4.25(l . By the resolvent identity, we have 

Ri K \z) - R € (z) = -R { e K) (z) (V^ - V) R e (z), R^\z) = (ffW - z)- 1 

with (V^ — V)R e (z) — > in operator norm. This shows that Ri K \z) — » R e {z) in operator norm. 
Using g^gt with k = for H ( , k) it is easy to check that H e Ri^ (i) is bounded in operator norm 
uniformly w.r.t. k. The resolvent idendity then shows that H e R^ K \i) — > H e R e (i) in operator 
norm and thus H t R i f'\z) converges as well. By induction H^ +1 R[ h '\z) k+1 — » H k+1 R e (z) k+1 in 
operator norm with a uniform bound of the form C(Re z) M /\lm z\ M (apply d\ to the resolvent 
identity for instance). Then, we see that {dlR { f' \z))R { ,? ] "(zy- 1 converges in S q /j with the same 
kind of bound as in lemma 14.111 and we can repeat the proof of proposition 14.121 The expected 
convergence follows from Helffer-Sjostrand formula and dominated convergence. □ 

5 The proof of theorem 11.21 

This section is entirely devoted to the proof of theorem 11.21 and more particularly to the proof 
of the asymptotic expansion. The continuity of £ q on the absolutely continuous spectrum is a 
consequence of the method. More precisely, we shall explicitely show that £ q (/i, h) is continuous in 
a neighborhood of fi = 1. The proof of the continuity of £ q (A), which we omit, would follow from 
the same method using (11.13(1 for fixed h. 

5.1 Isozaki-Kitada's method in the asymptotically hyperbolic case 

Let us consider a covering of Y by a finite number of open sets U n -\, each one of them being 
a relatively compact susbet of a coordinate patch U n -\ and consider their respective images on 



di(r) M (f(H e )-Q f Ne )(r) 
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M n_1 , i.e. £1 and f2, under the coordinates maps. For each such Q, we choose f2& for k = 1, • • • , 5 
such that 

Furthermore we can assume that Qk is open and convex (this will be usefull only for proposition 
15. 5|) for all k. Then we consider, the outgoing areas 1^ D D • ■ ■ D defined by 

r+ =T+(R k ,e k ,w k ,n k ), fc = l,-..,5. 

We assume that R > 1 and e,w G (0, f ) thus it is clear that C ri~_, for k = 2, • • • , 5. More 
precisely one can always find a cutoff function supported in Tt-i which is = f on Tt- We can 
choose for instance 

XB."(r)xn k (y)x+(p)xe"( e ~ 2r 9(y,v))xw"(p 2 + e~ 2r g{y,v)) (5-1) 

with smooth functions Xn k > Xn fc > X e fc an d Xu, fc such that 

supp XR k C (i? fe_1 ,oo) and X/f fc = 1 on (i? fc ,oo), 

supp xo fc C O fc _i and xn k = 1 on Sl t 

supp Xe fc C (— oo,e k ~ 1 ) and Xe fc = 1 on (~°°i £k )i 

supp Xii> fc C (1 - w k ~ 1 , 1 + w k ~ 1 ) and Xw k = 1 on (1 - w fe , 1 + w fc ). 

Since we will choose w; and £ small, it is enough to consider x+ supported into (0, oo) and such 
that x+ = 1 on (1/2, oo). Note that on Tt, we always have 

(l- w k -e k ) 1/2 <p<(l + w k ) 1/2 . (5.2) 

We are now ready to construct the functions needed for Isozaki-Kitada's method in the hyperbolic 
case. We start with the following proposition. 

Proposition 5.1. For R large enough, and e,w small enough, there exists a smooth function 
iff (r,y, p,rf) defined, for any e £ [0, 1], on such that 

(d r ipt) 2 + e- 2r g e (r,y,d y ipt) = p 2 . 

This function satisfies the following estimates for e £ [0, 1] and (r, y, p, i]) £ 

\did k d l p d^ { 9 t -rp- y.rj) | < C j<k<l , a<0 e-U +l >\ V \. (5.3) 

This proposition solves the equation l|2.22|l in the casep e (a;, £) = p 2 +e~ 2r g e (r, y, rj) and explains 
how differentiation w.r.t. e provides exponential decay. 

Proof. We follow the principle explained in subsection 12. 31 using the function S^(t, r, y, p, -q) given 
by proposition 13.61 Since solves H3.18|l and is a generating function of the flow (see (|2.2h'|l ) we 

have 

dtS+ = p e (r, y, d r St, d y St) = Pe(d P S?, d^S+, P , r,) 

= p 2 + e- 2d ^g t {d p St 1 d n St,r 1 ) 

where one remark that the last term of the second line is 0(e~ r ~ 2tp ) by Ij3.2t)|l hence is integrable. 
Thus we can use formula |2.28l with S+(t,p,rf) — tp 2 to define ipf, and then 1)5. 3J1 is a direct 
consequence of proposition 13.61 □ 

The next lemma is a preparation lemma for the resolution of the transport equations. We do 
not prove it since it can be obtained very similarly to the estimates on the geodesies of section [3] 
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Lemma 5.2. There exist R large enough and w,e small enough, the following result holds: for all 
(r, y, p, n) G and e G [0, 1], the solution (r*, y')(r, y, p, n) of 

2d r ^+(fl,yl,p,r 1 ) \ ff° £ \ /V 

is defined for all t > and satisfies (f * , y* , p, n) G . Furthermore, we have the estimates 

\dtd k r d l p dffi (ft-r- 2pt) | + \did k r d l p d^ (yt -y)\< C Ma ^-^ r \r,\ (5.4) 

for all e G [0, 1] and (r, y, p, n) G . 

This lemma, which gives in particular a precise sense to (|2.29|) in the current context, allows 
to define functions ai°\ • • • , ai of r, y, p, n G according to the formulas l|2.30|l . (|2.31|) where 
one has of course to replace x\ by f\ , y' and use the following explicit expression 

c e = d 2 r vt + e- 2r 9.(r,y,d y )vt + e- r £ ^(e^, v)^ dyY ^Wt > 

2+|a|=l 

where v'' Q = v 1 q U + t{v l { a — Wp"), with the notations of (|1.8|l . By l|5.3|) we see easily that 
c e (r,y, p,n) = 0(e~ r (i])) (note the factor (n) instead of \n\ which is caused by the term d l r ipf 
with I = 1 in the expression of c e ). This implies that the integral in Ij2.30|l is convergent since (15.211 
and l|5.4|) show that, for some c > 0, 

c £ (flyl,p,7 1 ) = 0(e~ ct - r {r 1 )), t>0. 

Thus ai is well defined on T^- By induction, one checks that, for m > 1, we have 

P £ (r, y, D r , D y )a^- 1] (f*. & P, l) = O^^in)), t > 0, 

on rj, for all e G [0, 1] and thus ai is well defined on for all m. More generally, by mean of 
proposition 15. II lemma l5~2l and (|5.2(l . we obtain easily the following result. 

Proposition 5.3. The functions ai°\ • • • , a\ N ^ defined in by ^2.3U$ and ^2.31\ satisfy 

d^dldffi (oi '(r, y, p, r?) - l) | < C ilWia ^e-C»' +1 ><» 7 > ) 

^9 r fc ^9«^a( m )(r,y,p,77)| < C^^e"^ 1 ^), l<m<N 

for all e G [0, 1] and (r, y, p, rj) G r \ . 

The details of the proof are left to the reader. They follow from the explicit expressions of the 
functions a[™'. 

In order to consider functions defined on R 2 ™, we choose a cutoff Xifik r i Vi Pi v) of the form (|5.1|l 
which is supported in and = 1 on . Then we multiply the functions ai°\ • • ■ , ai N ' by X2,3 
and define 

a e =X2,3 (aM+haP+-.- + h N aW). 
Then, following Isozaki-Kitada's method as explained in subsection !2.3l we consider 
H e J(<p+,a e ) - J((p+,a e )P, with P = h 2 D 2 r . 
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Of course, the choice of P follows from l|2.27[) since lim^oo p e (r, y, p,rj) = p 2 =: p(p). The operator 
H e J((p+, a e ) — J{<ft i a t)P has the following expression 

J(<p+,a! e ) + h N+2 j(<p+, X2 ,3Pe(r,y,D r ,D y )4 N ^ (5.5) 

where a' e is a linear combination of products of derivatives (of order > 1) of %2,3 and of derivatives 
of ai°^ , • • ■ , . The first term of l|5.5|l is produced by all the derivatives due to H e which may 

fall on X2.3- The amplitude of the second term is nothing but X2,s(o-e H h h N+2 ai N+2 ^ ) with 

the notations of subsection 12. 31 This follows from the construction of ipf and ■ • ■ , ai which 
satisfy respectively Hamilton- Jacobi's equation and the transport equations on the support of X2,3- 
The negligibility of i|5.5[) . or more precisely of (|2.17() . will be a consequence of the next lemma. 

Lemma 5.4. For R large enough, and e, w small enough we have the following property: for any 
symbol b(r, y, p, rf) supported in and such that, for all M , 

\d^d l p d^d^b(r,y,p, V )\ < C kMM (r)- M 

we have the following estimates for s > 

\\(r) M j(^,<)U( S )J^t,b)^r) A ^\ x <C M h M (s)- M , (5.6) 

\(r) M J^t,X2, 3 Pe(r,y,D r ,D y )4 N '>) U(s) J(<p+, 6)*<r) M | ^ < C M ^' (s)~ M (5.7) 

for all M and some universal constant no. Here U(s) — e^ % ^ p is the propagator of P. 

Note that the power h~ n ° on the right hand side of Ij5.7[l is harmless since we have a power 
h N+2 in H5.5J) . with N arbitrarily large. The proof of this lemma is not very hard and follows from 
suitable integrations by parts on the Schwartz kernels of the operators which are explicitly given 
by oscillatory integrals. However its proof is a bit long and we have postponed it to appendix 1X1 

The last step of the construction is the factorization of pseudo-differential operators. This is 
the purpose of the following proposition. 

Proposition 5.5. There exists R large enough and e,w small enough such that the following 
property holds: for any symbol c e € S, supported in T^f , one can find symbols bi°\ • • • , bi N ^ G S 
supported in such that 

J((pf,a e )J(tpf,b e )* — Oph(c e ) is h N negligible 

with b t = bf ] + hbi 1] + ... + h N bi N) . 

Proof. As explained in subsection l2.3l we need to study the map 

(p, 77) 1 ^ $+(r, y, r', y' , p,rf) = f d r , y ipt{r' + t{r - r'), y' + t(y - y'), p, rf) At (5.8) 



It is then clear that if {r,y,p,rf) and (r' ,y' , p,rf) belong to T + (R, e, w, Q) with f2 convex, then 
(r' + t(r — r'),y' + t(y — y'), p, rf) G T + (i?, Cos, w' , Vt) for all t e [0, 1], and any w' > w which can 
be chosen as close to w as we want by decreasing e. Here Co is given by Ij3.1|l . In particular, 



|3 7 ($+ - (fi,T))) \ < C* 7 e- (j+1)min(r ' r ' ) | 
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with d 1 = d{d^d^, d l p dydy, d@, for all (r, y, p, rj) and (r', y', p, rf) in T + (R, e, w, 0). Using the formula 
(|2.32(1 (with e = $+), we can define and by iterations b^\ • ■• , from c e . In particular, if c e is 
supported in T§, be,bi, ■ ■ ■ will be clearly supported in if R is large enough and e, w small 
enough. In order to check that the difference J(<pf ,a e )J(iff ,b € )* — Oph(c e ) is h N negligible, we 
look at its Schwartz kernel which we split into three terms using the partition of unit 

1 = 6 (r - r') + 6+{r - r') + 6+{r' - r) 

with 9q — 1 close to and 6 + supported in [1, 00). The term corresponding to #0 behaves nicely by 
construction. The off diagonal part is 0(/i°°) and is the sum of two oscillatory integrals which are 
either 0{e^^ +v ^ r ) or 0{e~^ +u ^ r ) after application d{, according to the sign of r — r' . We don't 
go any further into details. □ 



5.2 The dependence on k and the remainders 

Our next task is explain how to deal with 'remainders'. In subsection 12.11 or in section QJ we 
have seen why and how a lot of expansions of operators, in powers of h, could be obtained. It is 
now necessary to show why the contributions of the remainders of such expansions can indeed be 
neglected in the expected expansion. We introduce the notation =^y , which we will use extensively 
in the next subsection. Its meaning is the following. 



Definition 5.6. We write 



tr 



(if) (*,/*) 



=»i n 

— N u 



ifTe (t,h) is a family of trace class operators for n G (0, 1], whose trace is C q 1 with respect to 
e G [0, 1], measurable with respect to t el, and if there exists C/v iril , independent of h, such that 



lim (tr [T^ K) (t, h)) \ exists in S'(R t ) and belongs to L X (R, dt), 

kJ.0 IV /Jq 



hm{tr(^)(t,.))} q 



dt < C NiUl h 



N-ri! 



(5.9) 
(5.10) 



If the same result holds when integrating on [0, +00), then we will write =jy 1 instead of =jy 



If Se (t, h) is another family of trace class operators for K G (0, 1], then 

tr(TW (t,h)-SW(t,h)) 0, 



trftW (t,hj) tr (Sto{t,h) 

and similarly for =^+- 

Operators (or rather traces) satisfying 15.9(1 and 15.10fl are of interest since 



lim J- 1 



1 {tr (T £ W(a))} q = ^ 1 lim{tr (r e W(t,fc))} q = 0{h N - n ^) in C°(R). (5.11) 



In practice, it is enough to show, for instance, that lim K _ > o{7 1 i'^(t, /i)} q =: {Te"'(t, h)}q exists, 
in the weak topology of bounded operator and that 



,(0), 



{T^&h)}^ < C N h N ~ n ^{t), 
lim 1 tr{T e W(t,/i)} q = tr{7f)(i,/i)} q 



teR, he (0,1], 
teR, he (0,1], 



{T^(t,h)} c 



< Ch- ,l2 {t) n2 , t e R, h e (0,1], K g [0,1] 



(5.12) 
(5.13) 
(5.14) 
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for some L 1 function tp in (|5.12[l . and some 712 > in l|5.14ll . This last condition could be weakened 
since we could clearly allow C to depend arbitrarily on h. Usually, (|5.13|) and l|5.14|l are easy to 
check and the non trivial part of the job is to show l|5.12[) . 

Before giving explicit examples, let us explain how we will use definition 15.61 Recall that we 
are studying (|2.5(l which is nothing but 

e^*tt(f(H^)U^(t)V^f(H^ dt (5.15) 

with ff — f which we choose to be supported close to 1. The method is the following: if for any 
N large enough, we can find T^(t, h) such that 

tr(/(FW)C7 e W(t)yW/(ffW)) ee" 1 te(T<$(f,fc)) 

for some ri\ independent of N, and moreover such that {tr(T £ ( ^(i,/i))} q converges in <S'(R t ) as 
K i with a limit in L 1 (IR,dt) satisfying 

— J ^ e*"' ljm{ta (T$(t, ft)) } q dt = h~ n £ h*a k {n) + G(h N - n ) 

with Qfc G C°(7) and 0(h N ~ n ) understood in the topology of C°(J), then we get the existence of 
the expansion This is due to (|5.11(l since it implies that 

h) - -L J + ™ e *M* lim {tr (T$(t, h)) } ^ dt = inC° 

with N — n\ — 1 arbitrarily large if N is large. 

We shall now give explicit and useful examples of operators satisfying l|5.9|) and H5.10J) . To that 
end, we will use extensively the fact that for all / 6 Cq°(R), in the strong sense, 

U^(t)f(H^)^U e (t)f(He), KlQ 

with a norm bounded by sup |/|. We omit the proof of this easy fact since it follows by section 
0]or more directly, by Helffer-Sjostrand formula for instance. We will also use the following well 
known lemma, which is valid if B^ is a family of bounded operators and a family of S q . 

Lemma 5.7. If -> B^ strongly andT^ -> in S q then B^T^ -> B^T^ in S q . 

Let 1l[ K) N {h) = V (k) (j{H*) - Ql}^, with the notations of theorem Q3I where we include 
the depence on k. 

Proposition 5.8. Let I be an intervalle such that <|i._?gj) holds. Then all f,f S C^°(I), with 
ff = f, we have 

tr (f(HW)U( K \t)K { $ e {h)) =# 0, m = n + q. 
Proof. We only show (|5.12l) . We first rewrite the trace as 

tr ((r)-^f(H^)U^ (t)(r)- M (r) M TZ^ N (h)(r) M ) 
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and then apply df 1 . By Leibnitz rule and lemma 12.21 we obtain an explicit expression of this 
derivative. We note that, for any < k < q — 1, 

(r) M %- l - k KW(h)(r) M = 0(fc"-(»-*W*) in S q/(q _ k) (5.16) 

and that this estimate holds for k > 0, continuously w.r.t k and e, which will allow to let k J, 0. 
This can been seen with the same argument as the one used for the proof of lcmma l2~D in subsection 
14.31 On the other hand, using the notations of lemma [2^21 in which we now take k into account, 
we can write (r)~ M t/2.Wfy (r)~ M as the integral over F/ of 

h -j (v (K \r)- M d l °Si K)Jo (r) M ) ((r)- M UW' f °(t Q )(r)- M ) ({r) M V^ d 1 ^ (r) M ) • • • 

■ • • [{r)- M U^ J '{t)(r)- M ^j ((r} M d l e ^Si K ^+^{r)- M ^ . (5.17) 

Note that we used the notation 12. 81 Then, as before, it is not hard to check that 

(r) M {V^yd[S^{r) M {r) M V T d[S^<f(r) M e S q/(q _,_ T) 

for r = or 1. By lemma HT7I we can let k J, in l|5.17|l . and using i|1.5f) we get for k = 0, 



Ed (r) M dl- 1 - k 'R^ ) N (h) (r) M 1 1 1 < Ch N - n -jf[\\(r)- M Ui<(ti){T 



-Ml 



2=0 



Then the result follows from (|2.11|l . □ 
The same method can be applied to other kind of operators. We only explain which mod- 
ifications to do. For instance, if d'{lZ^ N {h) can be written in any chart as 1)4.16(1 x e~ r then 
everything works almost as well. This typically waht happens when one considers remainders of 
the pseudo-differential expansion of the adjoint of i/ K ) /(He^) or with negligible operator involved 
in proposition l5.5l Actually, in this case we don't have l|5.16ll anymore, but the exponential weights 
in l(4.16JI can be put on both sides of 1)5. 17J) which then becomes trace class and we can conclude 
similarly. This situation is also of interest when one has to consider the remainder involved in 
(12.14(1 . For the latter, the exponential weights are provided by the differentiation w.r.t. e in view 
of the explicit form of the symbols in term of the flow and theorem 13. 21 

We can also study the contribution of the right hand side of i|2.17[l using lemma lF^I Here there 
is one more integral over s but we clearly have L 1 estimates by lemma 15.41 and the convolution 
trick works again. Note finally that in this case we will choose n\ — n + q + tiq. 



5.3 The core of the proof 

Below, we will use extensively the notation of the previous subsection with m = n + q + uq , 
which is independent of N. 

We start from the semi-classical Fourier transform of (|2.5I) with / supported in [1 — w 5 /4, 1 + 
w 5 /4] and / = 1 close to 1. We furthermore assume that the coefficients of V (hence of V^) are 
all supported in {r > R}. We shall indicate how to modify the proof in the general case, however 
recall that it is the case for the principal symbol. 

Using theorem 14. 131 we get a pseudo-differential expansion of V^f(H { € K) ) whose symbols can 
be splited in two parts using a partition of unit associated to l|3.6Jl so that get, for any N, 

tr(/(ifW)C/i«)(f)yW/(ifW))=^ ]T h h tr(f(H^)U^(t)Q, h {^)) (5.18) 

k<N, n 
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where ^ € S^-°° and either supp7 e ( ^ C T+(i?, w 5 /2, a+, ft) or supp7 e ( ^ e T" (R, w 5 /2, er_, ft) 
for some f2 as described in the beginning of subsection l5.ll some crj- > and i? large enough. Let 
us consider only one term of the right hand side of (|5.18|) corresponding to a symbol supported in 
T + (R, w 5 /2, er + , f2) (the case — is similar). Then, using the trick l|2.12(l . lemma inl and Egorov's 
theorem we see that 

tr (f(H^)U^ (t)Op h ( 7 $)) =£ J] ^ tr (/( ff e (K) )^ (K) (5-19) 

z<jv 

with 7^' e Sl'~°° supported in T^. Here again, we study only one term of the right hand side of 
(I5.19|) . We split the integral corresponding to the inverse Fourier transform (|5.15|) into two parts 
and consider first 

— j e^Hv[f(H^)U^(t)Op h (if)) dt. (5.20) 

The interest of considering positive times and a symbol supported in is that we can use Q2.17H 
with the results of subsection l5.ll Using the notations of this subsection, we have 

tT(f(H^)U^(jk)Op H (^)) ^ + tr(/(ffW)J(^+-W,aW)C7(t)J^+'W,&W)*) (5.21) 

and, by centrality, we can rewrite the right hand side as 

tr[j^t' [K \b^rf(H^)J( V P {K \a^)U(t)) =% Yl h m tr(Op h (4%)U(t)) (5.22) 

m<N 

with Cefm € S\ ,0 ° . Now the contribution of each term of the right hand side of (|5.22|) to the 
expected asymptotic is easy to obtain since 

So 2^/ +0 ° etAI * { tr {°Ph{^%)U{t))}^t = ^l +0 °JJ f f c m (r,y,p,r ] )e^^ t drdydpdr ] dt 

thanks to the easy and crucial remark that 

c m =lim{ C W} 6S?'-». (5.23) 

The stationary phase theorem (in the variables t, p) yields easily the asymptotic of the last integral 
in integer powers of h with coefficients which are smooth functions of p in the neighborhood of 1. 
Note that we use the fact that p is close to 1 on the support of c m (see for instance Ij5.2|l l thus the 
Hessian of the phase (/i — p 2 )t is non degenerate w.r.t. t, p at the stationary point t = 0, p = /i 1 / 2 . 

We still have to explain how to deal with the contribution of negative times. Here we use the 
following trick due to Robert 

— J ei* tr [f{H^)U^{t)Op h {^)) dt = — J tr (f(H^)ui K \t)Op h (jif)*) dt 

which is a simple consequence of the fact that tr(A*) = tv(A) combined with the centrality of the 
trace. Then we can use the same method since 

tr (f(H^)U^(t)Op h (^fy) =£ J2 h m tr(f(H^)U^(t)Op h (^%)) 

rn<N 
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where the symbols je%, are derivatives of the complex conjugate of 7' fi . Since they are supported 
111 T+ as well, we can repeat the same method and the conclusion follows. 

This completes the proof under the condition that the coefficients of V vanish outside {r > R}. 
Otherwise, we proceed as follows. We write V = Vo + Voo with Voo supported in {r > R} and 
Vq compactly supported outside {r > R + 1}. The contribution of Vx ( or more precisely Vx , ) is 
treated as before and thus we only have to consider 

tr (f(H^)U^\t)V f(H^)) ■ 

Let K be the compact subset of T*X defined by the conditions r < R + 1 and \p c — 1 < u>s/4. 
By the non trapping condition, there exists T > and R' such that Hx4>i T (K) C {i? < r < R'}, 
if Hx is the projection onto the base. Then we choose ro such that TlxftiiK) is disjoint from 
{r > ro} for \t\ < T and this ensure the fact that = <j>\ for \t\ < T. Thus we can use the trick 
(|2"T!?|) with K e = Vb/(-ff e (K) ) combined with Egorov's theor em, and then repeat the same method. 
□ 



A Proof of lemma 15.41 

Let us first recall that for any smooth and bounded function a(x, x', £) defined on R 3 ™ the operator 
Ah, with Schwartz kernel Ah(x,x') defined as 



A h {x,x') = {2nh)- n / e^ x - x '^a(x,x',Z) d^ 



is bounded on L 2 (M. n ) and we have the following bound 



UhW^KC max sup d^,dla(x,x',0 , he (0,1] 



(A.l) 



for some constants C and no independent of a and h. This is the usual Calderon-Vaillancourt's 
theorem. This kind of operators preserve polynomial decay in the sense that, for any M € K we 
have 



x) M A h (x) 



-Ml 



< Cm max sup 

\a+a'+/3\<n M R3n 



d^d^a(x,x',0 , he (0,1] (A.2) 



for some Cm and um depending only on M. This is an easy consequence of IjA.ljl . 

Below, we shall use (jA.lll as follows. Assume that we have an operator Kt, with Schwartz 



kernel /CS of the form 



lC h (r,y,r',y') = (2irh)- n / / ei^'^+i^-y^^e^^' ^a s h (r,y,r' ,y' , p, V ) dpd V 



with a smooth function ip s real valued and such that, on the support of ai, 

Id 1 ^ 3 ^, y,r',y',p,r))\ < C 7 (s), 
for all 7. If we know moreover that for some N > and M > 0, we have 



\(Pa h (r,yy,y',p,v)\ < C 7 h N (s)- M (r)- M (r')- M 
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then we obtain the estimate 

1 1 (r) M K s h {r) M < Ch N - no {s)- M+no . (A.3) 

This follows simply from (|A.1|I by considering a — e^" < h al. 
For any a(r, y, p, rj), the kernel of J{ft 1 a )U (sjJift ■• b)* is 

{2iTh)- n J J e**'( r ' y ' r '*'- p - r >)a(r, y, p, 77)6(7-', y' , p, 77) dpdn (A.4) 

where the phase function $ £ is real valued and given by 

®t{r,y,r',y',p,r)) = <pf(r,y,p,i 1 )~sp 2 -Lp+(r',y',p,ri). 

We shall use extensively the fact that for any 7, d 1 (<p e (r, y, p, rj) — rp — y.rj) — 0(e k ) on T~£. This 
is a direct consequence of l|5.3|l and show that, if supp a C T^, supp b C , then 



d p $ s £ (r,y,r',y',p,r)) = r - r' - 2sp + 0(e 2 ) (A.5) 

d v $l(r,y,r',y', P ,rj) = y - y' + 0(e 2 ) (A.6) 

Recall moreover that, in view of 1)5.2(1 . there exists Ci, ci > such that 

ci<p<Ci, (A.7) 



on the support of a(r, y, p, 77)6(7"', t/', p, 77), if supp a C and supp C T±. 

Let us start the proof of (|5.6(l . Here we shall use the fact that a' e in r \ \ . More precisely, 
we need to study symbols a(r, y, p, rj) of the form 

{drX2,z{r, y, p, 77)) <9 7 a £ (r, y, p, 77), (eT^d™ X 2 ,3(r, y, p, 77)) <9 7 a £ (r, y, p, 77) 

with k > 1, |a| > 1. We first consider the symbols a involving 8^X2,3- 

1— //<9 r /aHs on Xi? 3 • 

In this case, we have R 2 < r < R 3 , hence we get the fast decay w.r.t. to r. Furthermore 
r' > R A on the support of b, thus we have r — r' < — R 4 + R 3 -C for R large, and we obtain 

d p $ s e (r,y,r',y',p,ri) < -1 - 2ps < -l-2 Cl s. (A.8) 

Thus we can integrate by part with h (9 p $J) _1 D p in (|A.4|I . and we get as many powers of h(s)^ 1 
as we want. 

2— If d r falls on x 6 3 - 

Then, we have e 3 < e~ 2r g(y, rj) < e 2 , thus we get 

e- 2r 'g(y\rj<e 4 ^e 3 <e' 2r g(y,r,). 
Since e~ 2r \t]\ 2 < Cqc~ 2t g{y',rj) and e~ 2r g(y,rj) < Coe~ 2r \r]\ 2 , this implies easily that 

e 2r ~ 2r ' < C 2 e. 

Thus r — r 1 = O(loge) <C 0. In particular r < r' and since we have as much powers of (r'} -1 as 
we want, we use the simple fact that (r)^')^ 1 is bounded to get as many powers of as we 

want. Furthermore, (|A.8|I still holds, for e small enough, thus we can integrate by part as before 
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and get as many powers of h(s}~ 1 as we wish. 
3— If d r falls on x w ' 3 ■ 

Here we have \p 2 + e~ 2r g(y, n) — 1| > w 3 . On the other hand we also have 



P 2 



+ e 2r g(y, V ) = p 2 + e 2r 'g{y', 77) + (eT^giy, 77) - e 2r ' g(y' , 77)) 

= p 2 + e- 2r 'g(y', V ) + 0(e 2 ) (A.9) 
with \p 2 + e~ 2r g(y', 77) — 1 1 < nr. Thus, on the support of a(r, y, p, rf)b(r', y', p, 77), we obtain 

\p 2 +e~ 2r g(y,ri) - 1| > w 3 and \p 2 + e- 2r g(y, 77) - 1| < w 4 + 0(e 2 ). (A.10) 



If we choose e < w 2 small enough, then IjA.lOJI can not hold and then a(r, y, p 7 rf)b(r', y 1 , p, 77) = 0. 
We continue our analysis by considering symbols involving e~' Q ' r 9^X2,3- 

4— If e~ r dy falls on \n 3 ■ 

Then y £ f^2 \ ^3 and y' € CI4 <s ^3, thus \y — y'\ > c > 0. Using (|A.6|I . with e small enough, 
we have > c/2 and we can integrate by part using h 2 \d v &l\~ 2 A v . This provides as many 

powers of h as we want. Furthermore the factor e~ r yields the fast decay w.r.t. r. We still need 
to explain how to get fast decay w.r.t. s. To that end, we introduce the following partition of unit 

-«-(^)+»°(^) +9 +(^) (A -"» 

Here 9- is supported in (—00, ci), 6 + in (3Ci, +00) and we can assume that, on the support of 6q, 
we have 

V — T 

ci/2< —— <4d. 
1 + s 

On this support, we can write 1 = (r — r') (r — r')~ M = (r — r') M (D((s)~ M ) , for any M. Since 
we already have fast decay w.r.t. r and r', the fast decay w.r.t. s follows. On the support of 
9-{{r - r')/(l + s)) (resp. 9+), using l|A.5|> and 1A.7I) , we see that, 

d p ^(r,y,r',y',p,r l )<-c 1 s + c 1 + 0(E 2 ) (resp. > C v s + 3d + 0{e 2 )) (A.12) 

thus, for s large enough, we can integrate by part as in 1— and get as many powers of (s) _1 as we 
want. 

5— If e~ r d y falls on \e 3 - We proceed as in 2— 

6— If e~ r d y falls on x w 3 ■ We proceed as in 3— 

This completes the proof of (|5.6() using (|A. 3|) . 

We now turn to the proof of Ij5.7|l . Here we just need to get as many powers of {r)^ 1 (s)^ 1 as 
we want. The estimates will rely upon the fact that 

d k r d l p d^dPp,{r, y,D r ,D v )aW\ < C^pe" (A.13) 

on for all k,l,a,/3, by proposition 15.31 We proceed as follows. We introduce the partition of 
unit (|A.lip and consider first what happens on the support oi 9o((r — r' ) / ( 1 + s) ) . Here we remark 
that 

e~ r (77) < e- ClS / 2 - r '(77) 
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where e r \rj\ is bounded on the support of b. This yields an exponential decay in s. Furthermore, 
we have r < r' + 4Cis, thus we can write for any M 

1 = (r)- M (r) M = (t)- m O (((r') + (s)) M ) . (A.14) 

The positive powers of (r') and (s) are respectively controlled by the fast decay w.r.t. r' of b and 
the exponential decay w.r.t. s, we obtain thus the expected decay. 

On the support of 9-{[r — t')/(1 + s)), we have similarly r — r' < Ci(l + s) and we can use 
(|A.14|I again. This yields (r)~ M for any M but we have to control (s) M . By the same method as 
in 4—, using the upper bound given by JUISJ), we get as many negative powers of (s) as we want. 

The last step is to study what happens on the support of 9+((r ~ r')/{\ + s)). We have 
r > r' + 3Ci(s + 1) and we get exponential decay in time, since 

e- r { V ) < e~ 3ClS ~ r '(77). 

Furthermore, by choosing e small enough, the lower bound in (|A.12fl shows that 

d p $ s e (r,y,r',y',p,ri)>Cis + 2C 1 >0. 

Thus we can integrate by part and get as many negative powers of d p <5>l as we wish. Then we note 
that 

\d p ^(r, y y,y',p,r,)\- M < C(r - r')- M (s) M < C'(r)- M (r') M (s) M 
and the fast decay with respect to r follows as before. This completes the proof. □ 

B Propagation estimates 

In this appendix, we give sufficient conditions leading to (|1.13|) . 

Proposition B.l. Assume that there are positive numbers w,ho,M,C such that 

sup \\(r)- M R e ( f x±iS)(r)- M \\ oo <Ch-\ (B.l) 

<5>0 

for all [i £ I, h e (0, ho] and e £ [0, 1] . Then holds. 

Proposition B.2. Assume that the manifold (X, G) is non trapping, and that the principal symbols 
of Pq and P\ coincide outside {r > ro}. Let I be a neighborhood of 1. Then there exists ro large 
enough, ho small enough and C > such that (If?, ifl holds with M = 1. 

The first proposition follows from Kato's theory of smooth perturbations. We recall its simple 
proof for the sake of completeness and to emphasize the uniformity with respect to the parameters. 
The proof of the second one uses Mourre theory and more particularly a combination of the 
ideasof [THj and HaEl. 

Proof of vrovosition WH Since we can always write / = /1/2 with /1, £ Cg°(J) and 

\\{r)- M f{H t )U e {t){r)- M \\oo < ||<r}- M / 1 (f/ £ )C/ e (i)||oc \\({r)- M f 2 {H £ )U € (-t)y\\ oc 

it is enough to estimate the norm of ||(r)" M /(iJ e )f7 £ (t)|| 00 in L 2 (R,dt). Let A e = (r)- M f(H e ), 
then by Parseval's identity (see we have, for all u s L 2 (X), 

2nh- 1 [ e- 25t ' h \\A e U € {t)u\\ 2 dt = [ \\A £ (R e (fi - iS) — R e {^ + iS))u\\ 2 dp,. (B.2) 
Jm Jr 
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On the other hand, (i/2)(R e (z) ~ R e (z)) — (Im z) R e (z)R e (z) defines a nonnegative operator for 
Im z > and we denote its positive square root by K e (z). We have trivially | \K e (z)A*u\ | 2 = 
(2iw)~ 1 (A*u, (Re(z) - R e (z))A* e u) and this yields 

ij \\AK £ (^ + i5) 2 u\\ 2 d^i^Chr 1 J \\K e (fj, + iS)u\\ 2 dfi = Ch- 1 7r\\u\\ 2 (B.3) 



since \\A e K e (/j, + iS) 2 u\\ < \\K e (fi + iS)A*\\\\K e (fj, + iS)u\\. Here C depends only on C in (|B.1|> and 
/, and is uniform w.r.t. e, h and S. The left hand side of (|B.3|) is nothing but the right hand side 
of (|B.2|) and the result follows easily. □ 

Proof of vrovosition \B.R The estimate (|B.1|) follows directly from the method of |37] provided the 
following Mourre estimate holds, with E e (I) the spectral projector of H e on /, 

E € (I)i [H e ,A h ] E e (I) > chE e (I) (B.4) 

for some constant c > independent of h e (0, ho] and e e [0, 1]. If (|B.4(1 holds then, the theory of 
Mourre shows that || \A h + i\~ 1 R e (n + iS)\A h + i^ 1 U = O^chy 1 ) thus (|B~T)) holds, provided 

||<r>- 1 A h || 00 < C u he(0,h o ]. (B.5) 

This reduces the proof to the construction of Ah- We now sketch its construction and refer to 
[TBI 1213 HHj for the details. The idea is to construct A h = A° + A™ with A^ supported near 
infinity and A^ compactly supported. Following JHj) w e define A^ as 

A< h = \f( H i) (0 2 oJ 2 s uhD r + hD r 6 2 uj 2 s u) /(Fx) (B.6) 

with / e C§°, f = 1 close to 1, and 9 = 9{r/R) smooth, bounded and supported near infinity, say 
in r > R/2. We also have 

u s = u( 2r - l0g{h * AY + 1) ), u = 2r + S-lo g (h 2 A Y + l), S > 0, 



where the function u> € C°°(K) is supported in [— 1, oo) and such that u> = 1 on [—1/2, oo). Then 
using the calculations of ^H] and pseudo-differential calculus, we see that 

f(H e )i[H e , AJ?]f(H e ) > h6u s f(H e )LJse + t(R, S)O(h) + 0(h 2 ) (B.7) 

if / is supported where / = 1. Here t(R, S) J. as R, S — > oo and the notation 0{h k ) holds in 
operator norm, uniformly w.r.t e. Note that the spectral cutoff f{H\) in l|B.6l) doesn't commute 
with H e and we use the fact, among other ones, that {f{H\) — f(H e ))9hD r = 0(R~°°). 

If there was no term u>s in the right hand side of (|B.7(I , we would have done half of our program. 
How to neglect 9(1 — ujs)f{H e )l We can give a pseudo-differential expansion of us similar to the 
one given in section^Jand thus, up to an operator which is 0(h), we can replace 1 — lus by a pseudo- 
differential operator with principal symbol 1 — u>((2r — log(g(y, rf) + l))/S). On its support, we have 
e~ 2r g(y,rj) + e~ 2r > e s ^ 2 . On the other hand, on the support of the principal symbol of f{H e ) 
we have e~ 2r g e (r,y,i]) < 3/2 (if / is supported close to 1) and thus e~ 2r g(y,-q) < 3/2 + 0(e~ r ). 
All this shows that the principal symbol of 9(1 — ojs)f(H e ) is identically for all R and S large 
enough, which implies that 

9(1 - uj s )f(h 2 P e ) - 0(h) (B.8) 
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where one should notice that O depends on R and S. Then we get 

f(H e )i[H e> A^]f{H e ) > h6f 2 (H £ )9 + l(R, S)O(h) + 0(h 2 ) (B.9) 

where 0(h 2 ) depends on R, S, but not O(h). 

We will now construct by the method of |2()l 14*5] , such that 

f(H £ )i[H e ,A° h }f(H e ) = h9f(H e )9 + Z(R,K)0(h) + 0(h 2 ), (B.10) 

with 9 E such that 8 2 + 9 2 = 1, and t(R,K) | as R, K f oo. Here K is another large 
parameter introduced below. If this holds, we easily get (|B.4|) by summing (|B.9|) and l|B.10|) . 
choosing first R, K and S large enough, then h small enough and then by multiplying both side 
by E e (I). Note that we also use the fact that [f(H e ),6] = 0(h). 

The idea is to define A^ as a (bounded) pseudo-differential operator whose principal symbol is 
the following function a, which is invariantly defined on T* X 

a =X / 2 o<p\dt fopi 
Jo 

where p\ and <j)\ are the principal symbol of H\ and the associated flow, and \ — x( r /K) is a 
Cq° function such that yj) 2 — 2 . Note that the integral is convergent thanks to the non trapping 
condition. Then the crucial remark is that the Poisson bracket is 

{ Pe , a} = f( Pl )8 2 + OiK- 1 ) + <D(R-°°) 

with 0(_fir _1 ) uniform w.r.t e (but not w.r.t R) in the topology of smooth and bounded functions. 
This follows first from the fact that 

{ Pl ,e 2 o 4>\} = { Pl ,e 2 } o # = | (e 2 o 4>\) 

and from the fact that p e — pi for r < tq with ro which we can choose > R. Then (|B.10J| follows 
from semi-classical pseudodifferential calculus, (see D 
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